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Abstract 

In our previous paper [12], a general framework was outhned to treat 
the approximate solutions of semilinear evolution equations; more precisely, 
a scheme was presented to infer from an approximate solution the existence 
(local or global in time) of an exact solution, and to estimate their distance. 
^ . In the first half of the present work the abstract framework of [12] is ex- 

! tended, so as to be applicable to evolutionary PDEs whose nonlinearities 

' contain derivatives in the space variables. In the second half of the paper this 

extended framework is applied to the incompressible Navier-Stokes equations, 
, on a torus T'^ of any dimension. In this way a number of results are obtained 

I in the setting of the Sobolev spaces EI"(T°'), choosing the approximate so- 

■ lutions in a number of different ways. With the simplest choices we recover 

local existence of the exact solution for arbitrary data and external forces, as 
well as global existence for small data and forces. With the supplementary 
^ ' assumption of exponential decay in time for the forces, the same decay law is 

^ . derived for the exact solution with small (zero mean) data and forces. The 

interval of existence for arbitrary data, the upper bounds on data and forces 
for global existence, and all estimates on the exponential decay of the exact 
solution are derived in a fully quantitative way (i.e., giving the values of all the 
necessary constants; this makes a difference with most of the previous litera- 
ture). Nextly, the Galerkin approximate solutions are considered and precise, 
still quantitative estimates are derived for their H" distance from the exact 
solution; these are global in time for small data and forces (with exponential 
time decay of the above distance, if the forces decay similarly). 
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1 Introduction. 



This is a continuation of our previous paper [12] on the approximate solutions of 
semihnear Cauchy problems in a Banach space F, and on their use to get fully 
quantitative estimates on: (i) the interval of existence of the exact solution; (ii) the 
distance at any time between the exact and the approximate solution. 
In [12], we mentioned the potential interest of (i) (ii) in relation to the equations 
of fluid dynamics. Here we treat speciflcally the incompressible Navier-Stokes (NS) 
equations on a torus T'^ of any dimension d ^ 2, taking for F a Sobolev space of 
vector fields over T*^. To be precise, we consider the Sobolev space EI"(T"') = H" 
of the "velocity fields" v : ^ 'R^ whose derivatives of order ^ n are square 
intcgrablc; then wc choose F := HI"^, where the subscripts so indicate the subspace 
of formed by the divergence free, zero mean velocity fields / : T"' ^ R'^ (of course, 
the condition of zero divergence represents incompressibility; the mean velocity can 
always be supposed to vanish, passing to a convenient moving frame). We always 
take n > d/2. 

The choice of T'* as a space domain allows a rather simple treatment, based on 
Fourier analysis; we presume that the results of this paper could be extended to 
bounded domains of R*^, with suitable boundary conditions. In our notations, the 
NS Cauchy problem is written 

cp{t) = A^{t) - £((^(t). d^{t)) + at) , m = /o , (1-1) 

where / = (p(t) is the velocity field at time t, £ the Leray projection on the diver- 
gence free vector fields and C,{t) is the external forcing at time t (more precisely, what 
remains of the external force field after applying £ and subtracting the mean value). 
We can regard (1.1) as a realization of the abstract semilinear Cauchy problem 

<f{t) = A^it) + V{^{t),t) , ^{to) = /o (1.2) 

where ^ : / i— > Af is a linear operator and V : (/, t) y-^ V{f, t) is a nonlinear map. 
Of course, in the NS case A is the Laplacian A and V{f,t) := —£,{f»df) + ^1). By 
a standard method, both (1.1) and its abstract version (1.2) can be reformulated as 
a Volterra integral equation, involving the semigroup (e*'^)t^o- 
In [12], a general setting was proposed for semilinear Volterra problems, when the 
nonlinearity V{-,t) is a sufficiently smooth map of a Banach space F into itself. This 
setting cannot be applied to Cauchy problems like (1.1). In fact, due to the presence 
of the derivatives df, the map / t— > £(/• df ) cannot be seen as a smooth map of a 
Sobolev space, say H^g, into itself; on the contrary, the above map is smooth from 
to IHI^q"^. The external forcing ^ : t ^-^ ^t) fits well to this situation if we require 
it to be a sufficiently smooth map from [0, +00) to H""-'^. 

In view of the applications to (1.1), in the first half of the present paper (Sections 
2-5) we extend the abstract framework of [12] to the case where, at each time t, 
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V{-, t) is a smooth map between F and a larger Banach space F_. A general scheme 
to treat approximate solutions is developed along these lines; this could be applied 

not only to (1.1), but also to other evolutionary PDEs (essentially, of parabolic type) 
with space derivatives in the nonlinear part. 

In the second half of the paper (Sections 6-10) we fix the attention on the NS 
equations, in the framework of the above mentioned spaces (incidentally, we 
wish to point out that other function spaces could be used to analyse the same 
equations within our general scheme). 

Some technicalities related to either the first or the second half are presented in 
Appendices A-H. 

Of course, there is an enormous literature on NS equations, their approximation 
methods and the intervals of existence of the exact solutions: references [2] [3] [5] [6] 
[7] [8] [9] [10] [11] [14] [15] are examples including seminal works, classical treatises 
and recent contributions. Some differences between the present analysis and most 
of the published literature are the following: 

(i) Our discussion of the NS approximate solutions is part of a more general frame- 
work, in the spirit of the first half of the paper. 

(ii) Our analysis is fully quantitative: any function, numerical constant, etc., ap- 
pearing in our estimates on the solutions is given explicitly. In the end, our approach 
gives bounds on the interval of existence of the exact NS solution and on its dis- 
tance from the approximate solution in terms of fully computable numbers; such 
computations are exemplified in a number of cases. 

(iii) If compared with other contributions, our approach seems to be more suitable 
to derive the existence of global exact solutions from suitable approximate solutions, 
under specific conditions (typically, of small initial data); a comment on this point 
appears in Remark 8.7 (iii). 

Hereafter we give more details about the contents of the paper. 

First half: a general setting for the approximate solutions of (1.2). We have 
just mentioned the assumption P(-, t) : F i— > F_. We further ly suppose ^ : F+ — > F_ 
where F_,_ is a dense subspace of F, to be equipped with the graph norm of A; in 
the end, this gives a triple of spaces F_,_ C F C F_. 

To go on, we require A to generate a semigroup on F_, with the fundamental regu- 
larizing property e*'^(F_) C F for all t > 0. A more precise description of all these 
assumptions is given in Section 2: here we suppose, amongst else, the availability of 
an upper bound u_{t) e (0, -|-oo) for the operator norm of e*-^, regarding the latter 
as a map from F_ to F. The bound u_ is allowed to diverge (mildly) for t ^ 0^, an 
indication that e*'^F_ (t F for t = 0: the precise assumption is u^(t) = 0{l/t^~'^) 
with < cr ^ 1. In apphcations to the NS system, F = and = H^^^; 
the semigroup (e*"^) of the Laplacian has the prescribed regularizing features, with 
a = 1/2. 
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In Section 3 we present a general theory of the approximate solutions, for an abstract 
Cauchy (or Volterra) problem of the type sketched above. The basic idea is to 

associate to any approximate solution ip^^ : [to,T) — F of the problem an integral 
control inequality for an unknown function TZ : [to, T) [0, +00); this has the form 



where £ : [to^T) — > [0, +00) is an estimator for the (integral) error of 99 and £ is 
a function describing the growth of V from (/7„p. The main result in this framework 
is the following: if the control inequality is fulfilled by some function TZ on \tQ,T), 
then the semilinear Volterra problem has an exact solution (/? : [to, T") — > F, and 



for all t in this interval (|| || is the norm of F). 

When F_ = F, we recover from here the framework of [12]. Similarly to the result 
of [12], the present theorem about 7^, (pap and (p can be considered as the abstract 
and unifying form of many statements, appearing in the literature about specific 
systems. 

The available literature would suggest to prove the above theorem along this path: 
(i) derive an existence theorem for (p on small intervals; (ii) use some nonlinear 
Gronwall lemma to prove that ||(^(i)— (^„j,(t)|| ^ 7^(t) on any interval [to, T') C [to, T) 
where </? is defined; (iii) show the existence of </p on the full domain [to, 7") of by 
the following reductio ad absurduni: if not so, Wfit) — </7„p(t)|| would diverge before 
T and its upper bound via 7^(t) would be violated. 

Our proof of the theorem on TZ, (p^^ and (p, presented in Section 4, replaces the 
above strategy with a more constructive approach. The main idea is to interpret 
the control inequality (1.3) as individuating a tube of radius TZ — TZ{t) around (^„j„ 
invariant under the action of the semilinear Volterra operator J for our problem. 
This makes possible to construct the solution by an iteration of Peano-Picard type, 
starting from p^^, the result is a Cauchy sequence of functions p^ = J^iVap) on 
[to,T), {k = 0,1,2,...), whose k — > +00 limit is an exact solution of the given 
Volterra problem. 

Prom this viewpoint, existence of the solution on a short time interval, with any 
datum p(to) = /o, is a very simple corollary of the previous theorem based the 
choice Pap(t) := constant = /o. 

Even though there is a basic analogy with [12], proving the main theorem on ap- 
proximate solutions is technically more difficult in the present case, mainly due to 
the divergence of 'U-(t) for t 0+. Such a divergence is also relevant in applications: 
in fact, differently from [12], Eq. (1.3) with ^ replaced by = cannot be reduced to 
an ordinary differential equation. Our assumption u_{t) = 0(l/t^~'^) relates (1.3) 




to 



(1.3) 



\\p{t)-pum^nt) 
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to the framework of singular integral equations of fractional type (which could be 
interpreted in terms of the so-called "fractional differential calculus"). 

In spite of these pathologies, solving (1.3) is rather simple when the semigroup 
(e*-^) and ip^j, have suitable features, and the nonlinear function V has the (affine) 
quadratic structure 

v{f,t) = nfj)+m , (1-4) 

with T : F_ X F_ — > F a continuous bilinear form and ^ : [0, -\-oo) — > F_ a (locally 
Lipschitz) map; this is the subject of Section 5 (where the datum /o of (1.2) is always 
specified at tQ — Qi). 

The section starts from a fairly general statement on the control inequaility (1.3), 
which is subsequently applied with specific choices of the approximate solution. 
First of all, we consider the choice (Pap{t) '■= 0. In this case, for any datum /o and 
external forcing ^, we construct for the control inequality a solution TZ with domain a 
suitable interval [0, T); this imphes the existence on [0, T) of the solution (p of (1.2), 
and gives an estimate ||v'(t)|| ^ T^{t) on the same interval. If /o and ^ axe sufficiently 
small, T = +00 and so (/? is global. With the stronger assumption that ^{t) decays 
exponentially for t — > +oo, we derive for the control inequality a solution t t— > TZ{t) 
which is also exponentially decaying; so, the same can be said for ||(/9(t)||. Next 
we consider, for a small /o and a small, exponentially decaying ^, the approximate 
solution (f^p obtained solving the linear Cauchy problem ^ap{t) = AiPap{t) + C,{t), 
V'ap(O) = /o- In this case the control equation still possesses a global, exponentially 
decaying solution TZ, giving a precise estimate on the distance ||v^(t) — ip^p(t)\\. 

Second half: applications to the NS equations. In Section 6 we review the 
Sobolev spaces of vector fields on T*^, and the Leray formulation of the incompress- 
ible NS equations within this framework; furthermore, we show that the Cauchy 
problem with mean initial velocity mo can be reduced to an equivalent Cauchy prob- 
lem where the initial velocity has zero mean, by a change of space-time coordinates 
{x,t) 1-^ (,T — h{t),t), where the function t i-^ h{t) is suitably determined. 
In the same section we give explicitly the constants Knd = such that H/.Ss'll^-i ^ 
-f^n||/||n||fl'||n for all vclocity fields f,g on T'*, || ||„ and || ||n_i denoting the Sobolev 
norms of orders n and n — 1. The study of these constants, inspired by our previous 
work [13], prepares the fully quantitative application of the methods presented in 
the first half of the paper. 

Section 7 starts from the formulation (1.1) of the Cauchy problem, in the already 
mentioned Sobolev spaces F = W^^, F^ = H^^\ We check that (1.1) fulfills all 
requirements of the general theory for quadratic nonlinearities, and construct the 

estimator u_ for the semigroup (e*"^). 

In Section 8 we rephrase for the NS equations all the results of Section 5 on the 
abstract quadratic case (1.4). The estimates on the time of existence T, for arbitrary 
data and forcing, have a fully explicit form; the same happens for the bounds on 
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the norms ||/o||n, ||C(^)||n-i which ensure global existence and, possibly, exponential 
decay of (p{t) for t — > +00. 

In Section 9 we discuss the approximate NS solutions provided by the Galerkin 
method. More precisely, for each finite set 0) of wave vectors we consider the 
subspace spanned by the exponentials e*^'"* {k E G), and the projection on of 
the NS Cauchy problem; this has a solution 1 1— > (p'^{t) (in general, on a sufficiently 
small interval; with special assumptions, also involving the forcing, (p*^ is global 
and decays exponentially for t — > +cxd). Applying the framework of Section 7 with 
<^ap = we derive the following results (with p > n and |G| := inf^g2;'*\G a/I + l^p)- 

(i) For any initial datum /q G M^q of the NS Cauchy problem (1.1), and each external 
forcing ^ with values in H^o ^, both ip'^ and (p exist on a suitable interval [0, T), and 
there is an estimate 

Mt)-^^WIIn^^^^f|^ fortG[0,T); 

T can be +00, if the datum and the forcing are sufficiently small. Both T and the 

function t ^-^ yVnp\G\{t) ^'I'c given explicitly 

(ii) If /o G is sufficiently small and there is a small, exponentially decaying 
forcing t ^ ^(t) G H|o \ then 

Mt) - ^^{t)\\n ^ l^e-* for t G [0,+oo) ; 

the upper bounds for /q, ^ and the coefficient IVfjp|G| are also given explicitly. 
The results (i) (ii) imply convergence of p'-^ to p as |G| +00, on the time interval 
where the previous estimates hold (which can be [0, +00), as pointed out). 
In Section 10 we exemplify our estimates giving the numerical values of T and of 
the error estimators in (i) (ii) for certain data and forcing, with d — 3 and n — 2, 
p = 4 . 



2 Introducing the abstract setting. 

Notations, (i) All Banach spaces considered in this paper are over the same field, 
which can be R or C. 

(ii) If X and Y are Banach spaces, we write 

X-^Y (2.1) 

to indicate that X is a dense vector subspace of Y and that its natural inclusion 
into Y is continuous (i.e., ||a;||Y ^ constant \\x\\^ for aU x G X). 

(iii) Consider two sets 9, X and a function x : 9 — > X, 1 1— > x(i). The graph of x is 

g^X-={ix{t),t) |tG9}cXx9. (2.2) 
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If X = [0, +00] , we define the subgraph of x as 

sgrx := {(r, t) | t G 9, r G [0, x{t))} C [0, +00) x 9 . (2.3) 

(iv) Consider a function x : 9 — X, where 9 is a real interval and X a Banach 
space. This function is locally Lipschitz if, for each compact subset / of 9, there is 
a constant M = M{I) e [0, +00) such that 

\\x(t)-x{t')\\^^M\t-t'\ ior al\t,t' e I . (2.4) 

As usually, we denote with C°'^(9,X) the set of these functions. 
General assumptions. Throughout the section, we will consider a set 

iF^,F,F_,A,u,u_,V) (2.5) 

with the following properties. 

(PI) F.,., F and F_ are Banach spaces with norms || || + , || || and || ||_, such that 

F+ F F_ . (2.6) 

Here and in the sequel, B(/o,r) will denote the open ball {/ e F | ||/ — /o|| < r} 
(the radius r can be +00, and in this case B(/o, r) = F). 
(P2) ^ is a linear operator such that 

^ : F+ ^ F_ , f^Af . (2.7) 

Viewing F+ as a subspace of F_, the norm || ||+ is equivalent to the graph norm 

/eF,^||/||_ + p/||_. 

(P3) Viewing ^ as a densely defined linear operator in F_, it is assumed that A 
generates a strongly continuous semigroup (e*-^)tg[o,+oo) on F_ (of course, from the 
standard theory of linear semigroups, we have e*'^(F_,.) C F.,. for all t ^ 0). 

(P4) One has 

e*^(F) C F for t e [0, +00) ; (2.8) 

the function (/, t) ^ e*"^/ gives a strongly continuous semigroup on F (i.e., it is 
continuous from F x [0, +00) to F). 

Furthermore, ueC([0,+oo),(0,+oo))isa function such that 

\\e'^f\\<u{t)\\f\\ fort^O, /eF; (2.9) 

this function will be referred to as an estimator for the semigroup (e*-^) with respect 
to the norm of F. 

(P5) One has 

e*-^(F_) C F for t e (0, +00) ; (2.10) 



7 



the function (/, t) i— > e*-^/ is continuous from F_ x (0, +00) to F (in a few words: 
for all t > 0, e*"^ regularizes the vectors of F_, sending them into F continuously). 
Furthermore, u_ e C((0, +00), (0, +00)) is a function such that 

||e*-^/|K«_(t)||/||_ fori>0,/eF_ ; (2.11) 

u_{t)^0{:^) forf^O+, ae(0,l]. (2.12) 

The function u_ will be referred to as an estimator for the semigroup e*"^ with 
respect to the norms of F and F_; Eq. (2.12) ensures its integrability in any right 
neighbourhood oit — 0. 

(P6) One has 

: L'omT' CF X R^F_ , {f,t)^V{f,t), (2.13) 

and the domain of V is semi-open in F x R: by this we mean that, for any (/q, to) G 
DomV, there are 5, r G (0, +00] such that B(/o, r) x [to, to + 5) C DomV. Further- 
more, V is Lipschitz on each closed, bounded subset C of F x R such that C C DomV; 
by this, we mean that there are constants L — L{C) and M — M{C) G [0, +00) such 
that 

\\r(f, t) - V(f', t') ||_ ^ L\\f - /'II +M\t- t'\ for all (/, t), (/', t') E C . (2.14) 

2.1 Remark. As anticipated, our aim is to discuss the Cauchy problem ip{t) — 
A(p{t) + V{(p{t),t), (p{to) = /o (and its equivalent formulation as a Volterra prob- 
lem) for a system (F+, F, F_,A, u, u_, V) with the previously mentioned properties 
(P1),...,(P6). In comparison with the present work, the analysis of [12] corresponds 
to the special case 

F_ = F, u_=u (2.15) 
in which, by the continuity of m at t = 0, Eq. (2.12) is fulfilled with cr = 1 {}). 

Preliminciries to the analysis of the Cauchy and Volterra problems. 

(i) In the sequel, whenever we consider an interval \tQ,T), we intend —00 < < 

(2). 

(ii) Let us consider a function uj G C{[to,T),F_) and the function 

n-.te [to,T) ^ n{t) := f dse^^-'^^u;(s) . (2.16) 

J to 

■^In [12] e*-^ was written U{t), and F+ was simply indicated with DomA; furthermore, we 
assumed DomP to be open in R x F. 

^In [12], we also considered solutions of the Cauchy or Volterra problems with domain a closed, 
bounded interval [to, T]; the symbol [to, T\ was employed to denote an interval of either type. Here 
we only consider the first case (semiopen, possibly unbounded), simply to avoid tedious distinctions. 
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Using the regularising properties (P5) of the semigroup with respect to the spaces F 
and F_, one easily proves that for each fixed t the function s i— > e^*~^^-^uj{s) belongs 
to L^{{to,t),dt,F) and n e C{[to,T),F). 

(iii) All the results in (ii) apply in particular to the case uj{s) := V{ip{s),s) where 
ip e C{[to, T), F) and gr^^ C DomV. Functions of this form will often appear in the 
forthcoming analysis of Volterra problems. 

(iv) Many facts stated in the sequel depend on the basic identity, here recalled for 
future citation, 

= e('-*°)-^V'(to) + f ds e^*-^)-^ U(s) - Ail^{s)\ , (2.17) 

holding for any function G C([to, T), F+) fl C^([to, T), F_) and each t in this 
interval. 

Formal definitions of the Cauchy and Volterra problems. These definitions 
are similar to the ones adopted in [12], with shght changes due to the present use of 
two different spaces F, F_. 

2.2 Definition. Consider a pair {fo^h) e DomV, with fo e F+. The Cauchy 
problem CP(/o,io) with datum fo at time to is the following one: 

Find if e C{[to, T), F+) n C^([to, T), F„) such that gnp C DomV and 

ip{t) = A^{t) + V{ip{t),t) for all t e [to, T) , ip{to) = fo ■ (2.18) 

Wc note that C([to, T), F+) C C([to, T), F). This fact, with the properties of A and 
V, implies the following: if (/? G C([to, 7"), F+) and gap C DomV, the right hand 
side of the differential equation in (2.18) defines a function in C([to, T), F_). 

2.3 Definition. Consider a pair [fo, to) & DomV . The Volterra problem WP{fo, to) 
with datum fo at time to is the following one: 

Find (f e C{[to,T),F) such that gnp C DomV and 
ip(t) = e(*-*°)-^/o + [ ds e^^-'^^V{ip{s), s) for aU t e [to, T) . (2.19) 

2.4 Proposition. For (/o,to) G DomV and fo G -F+, we have the following. 

(i) a solution p of CP{fo,to) is also solution o/VP(/o,to) ; 

(ii) a solution p o/VP(/o,to) is also a solution o/CP(/o,to) ; if F_ is reflexive. 

Proof It is based on (2.17): see [1]. The derivation of (ii), which is the most 
technical part, uses the Lipschitz property (P6) of V and the reflexivity of F_ to 
show that a solution of VP(/o, to) has the necessary regularity to fulfill CP(/o, to) ■ □ 
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2.5 Proposition. (Uniqueness theorem for the Volterra problem). Consider a 
pair (/o, to) ^ DomV , and assume that VP(/o, to) has two solutions (p e C{[to,T), F), 
if' eC{[to,T'),F). Then 

ip{t) = ip'{t) for t e [to, min(r, T')) . (2.20) 

Proof. We consider any r G [to, niin(T, T')), and show that = in [to,T]. To 

this purpose, we subtract Eq.(2.19) for if from the analogous equation for (/?'; taking 
the norm || || and using Eqs. (2.11) (2.12) (2.14), for each t e [to,T] we obtain: 

Mt)-v'{t)\\ ^ Jdsu_(t - s)\\P(^{s),s) -V(^'{s),s)\\_ (2.21) 
<UL 



In the above: L ^ is a constant fulfiUing the Lipschitz condition (2.14) for V 
on the set C := gT{ip \ [to,T]) U gr(9?' [ [to,r]); ^ is a constant such that 
u_{t') ^ U/t'^~'' for all t' e (0,t] (which exists due to (2.12)). 
Eq. (2.21) implies \\<f{t) — ^p'{t)\\ = for all t e [to,T]; in fact, this result follows 
applying to the function z{t) :— \\(f{t) — (f'{t)\\ the forthcoming Lemma. □ 



2.6 Lemma. Consider a function z G C([to,T], [O, +oo)) (with — oo < to < r < 
+oo), and assume there are A e [0, +oo), a e (0, 1] such that 

z{t) f ds for t e [to, r] . (2.22) 

J to l^- ■SJ 

Then, z(t) = for all t e [to,7"]- 

Proof. See Appendix A. □ 



2.7 Remcirk. For VP(/o,to)we will grant existence as well, on sufficiently small 
time intervals (see the forthcoming Proposition 3.10, where local existence is ob- 
tained as a simple application of the general theory of approximate solutions) . 

The Volterra integral operator. This is the (nonhnear) integral operator ap- 
pearing in problem VP(/o,to) . More precisely, let us state the following. 

2.8 Definition. Let (/o,to) G DomV. The Volterra integral operator Jt^f^^to) = J 
associated to this pair is the following map: 

(i) Dam J' is made of the functions ip e C{[to,T), F) (with arbitrary T e {to, +oo]) 
such that gnp C DomV; 

(a) for each ip in this domain, J{^) G C{\to,T), F) is the function 

t e [to, T) ^ Jim) e(*-*°)'^/o + /* ds e^'-^^^Vim, s) ■ (2.23) 

J to 
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3 Approximate solutions of the Volterra and Cau- 
chy problems: the main result. 

Throughout the section, we consider again a set {F^,F,F_, A,u,u^,V), with the 
properties (P1)-(P6) of the previous section. The definitions that follow generalize 
similar notions, introduced in [12]. 

Approximate solutions, and their errors. We introduce them in the following 
way. 

3.1 Definition. Let (/o,to) £ DomV. 

(i) An approximate solution of VP{fQ,to) is any function ip^p G C{\tQ,T), F) such 
that gripap C DomV. 

(a) The integral error of (pap is the function 

E{^J:=ip^,-J{ipJeC{[to,T),F} , (3.1) 

z.e., E{^J{t) = - e(*-*")-^/o - ^ ds e(*-)-^P(v.„,(s), s). 

An integral error estimator for (p^^ is a function E e C([to, T), [0, +oo)) such that, 

for all t in this interval, 

\\E{cpJ{t)\\^£{t) . (3.2) 

3.2 Definition. Let (/o,to) ^ DomV, and fo G -F+. 

(i) An approximate solution of CP{fo,to) is any function ip^^ G C{[to,T), F_^.) fl 
C^{[to,T),F_) such that gnp,^ C DomV. 

(a) The datum error for </?„p is the difference 

diipj := ifM - /o e C F ; (3.3) 

a datum error estimator for ip^p is a nonnegative real number 6 such that 

MvJW^S. (3.4) 

(Hi) The differential error of (p^p is the function 

e{^J e C{[to,T), F_), t ^ e{^J{t) := ^^.(t) - A^^,{t) - V{^^,{t),t) ; (3.5) 

a differential error estimator for tp^p is a function e G C{[to, T), [0, +00)) such that, 
for t in this interval, 

\\e{<pj{t)\\_^e{t) . (3.6) 

3.3 Remarks, (i) A function ip^p as in Definition 3.1 (resp.. Definition 3.2) is a 
solution of VP(/o, to) (resp., of CP(/o, to) ) if and only if E{ip^p) = (resp., d{(pap) = 
and e{(p^p) = 0). 

(ii) Of course, the previous definitions of the error estimators can be fulfilled setting 
£{t) ■■= m^apmi S := Mv^Jl e{t) := ||e((^.J(t)||_. 



11 



3.4 Lemma. Let {fo,to) e DomV, fo e F_^., and ip^p be an approximate solution 
of CP {fo, to) with datum and differential errors d{ipap), ^{fap}- Then: 

(i) fap is also an approximate solution o/VP(/o,to) ; with integral error 

E{^J{t) = e(*-*°)-^ d{^J + f ds e(*-^)-^e((^„J(s) . (3.7) 

(a) If S,e are datum and differential error estimators for ipap, an integral error 
estimator for (p^p is 

E{t) := u{t -to) 5+ f ds u_{t - s)e{s) for all t e [tQ,T). (3.8) 

Proof, (i) To derive Eq.(3.7), use the definitions of E{ip^p), d{ipap), e^ip^j) and tlie 
identity (2.17) with := ipap- 

(ii) To derive the estimator (3.8), apply || || to both sides of (3.7), using Eqs. (2.9), 
(2.11) for u, u_, and Eqs. (3.4), (3.6) for S, e. □ 

Growth of 7-* from a curve. To introduce this notion, we need some notations. 
Let us consider a function p e C([to, T), (0, +oo]); we recall that, according to (2.3), 
the subgraph of p is the set sgrp := {{r,t) \ t E [to,T) , r E [0, p{t))}. Furthermore, 
let (f) e C{[to, T"), F). We define the p-tube around cf) as the set 

T(0,p) := {{f,t) \feF,te [to,T), 11/ - m\ < Pit)} ; (3-9) 

of course, the above tube is the whole space F if p{t) = +oo for all t (^). 

3.5 Definition. Let </> G C{[tQ,T), F), with grcf) C DomV. A growth estimator 

for V from (j) is a function £ with these features. 

(i) The domain of £ is the subgraph of some function p G C([to, T), (0, +oo]), and 

i e C{sgrp, [0, +oo)) , (r, t) ^ £(r, t) ; (3.10) 

£ is nondecreasing in the first variable, i.e., £{r, t) ^ £{r', t) for r ^ r' and any t. 

(ii) The function p in (i) is such that T(0, p) C DomV. For all {f,t) G T(0, p), it 
is 

\\nf,t)-vm,t)\\-<mf-m\,t) . (3.11) 

3.6 Remark. Consider any tube T((;A, p') C DomV. Using the Lipschitz property 
(2.14) of V, one can easily construct a growth estimator £ of domain sgr(p'/2), 

depending linearly on r: ('(r.f) = X(t)r. 

[12], this notion was presented in the case p = constant; the present generaHzation is 
harmless, and could have been employed in our previous work as well. 
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The main result on approximate solutions. This is contained in the following 

3.7 Proposition. Let{fo,to) eDomV, and consider the problem \fP{fo, to) . Sup- 
pose that: 

(i) ifap e C{[tQ,T), F) is an approximate solution o/VP(/o,to)7 ^ £ C{[to,T), 
[0, +00)) is an estimator for the integral error E{(p^p); 

ii) i G C{sgrp, [0,+oo)) is a growth estimator for V from (p^j, (for a suitable p e 

C([to,T),(0,+oo]). 

Consider the following problem: 

Find n e C([to, T), [0, +00)) such that grTZ C sgrp, and 

E{t) + [ ds u_{t - s) £{TZ{s), s) ^ 7^(^) for t e [to, T) . (3.12) 

JtQ 

If (3.12) has a solution TZ on [to,T), then VP(/o,to) has a solution tp with the same 
domain, and 

Mt) - <PaM\ ^ nt) forte [to,T) . (3.13) 
The solution (p is constructed by a Peano-Picard iteration of J , starting from </7„p. 
Proof. See the next section. □ 



3.8 Definition. Eq.(3.12) will be referred to as the control inequahty. 

3.9 Remarks, (i) It is worthwhile stressing the following: the estimators u_,i,S 
in the control inequality (3.12) depend on A,V,Pap, and should be regarded as 
known when the Volterra problem and the approximate solution are specified. So, 
(3.12) is a problem in one unknown TZ, that one sets up using only informations 
about p^p. After TZ has been found, it is possible to draw conclusions about the 
(exact) solution p of VP(/o, to) • In the usual language: the control inequality allows 
predictions on p through an a posteriori analysis of p^p- 

(ii) (extending to the present framework a comment in [12]). Typically, one meets 
this situation: p^^p, £, £ are defined for t in some interval [to,T'), and the control 
inequality (3.12) has a solution TZ on an interval [to,T) C [to,T'); in this case one 
renames p^p, S, etc. the restrictions of the previous functions to [to,T), and applies 
Proposition 3.7 to them (as an example, this occurs essentially in the proof of the 
forthcoming result). 

A first implication of Proposition 3.7: local existence. The most general 
and simple consequence of Proposition 3.7 is the fact anticipated in Remark 2.7, 
i.e., the local existence for the Volterra problem. Here we formulate this statement 
precisely. 
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3.10 Proposition. Let (/o,^o) £ DomV. Then, there are R',T',£,£ such that 

(i)-(iii) hold: 

(i) R' e (0, +00], r e {to, +00] and B{fo, B!) x [to,T') C DomV; 

(ii) 8 e C([to, T'), [0, +00)) and, for all t e [to, T'), 

ll/o - e(*-*°)-^/o - f ds e^'-'^Mfo,s) \\ ^ S{t) , £{to) = ; (3.14) 

J to 

(Hi) £ G C([0, R') X [to, T'), [0, +cxd)), (r, t) 1— > £(r, t); t/iis function is non decreasing 
in the first variable and, for (/, t) e -B(/o, -R') x \to, T'), 

\\vu,t) -nhM- - hit) . (3.15) 

Given R',T',S,£ with properties (i)-(iii), we have (a)(b): 

(a) there are R e (0, R') and T e {to, T'] such that, for all t e [to, T), 

S{t) + / dsu_{t- s)£{R, s) ^R; (3.16) 

J to 

(b) if T and R are as in item (a), VP(/o,to) has a solution ip of domain [to,T) and, 
for all t in this interval, 

y{t)-fo\\^R. (3.17) 

Proof. Step 1. Existence of R',T',S,£. A pair {R',T') as in (i) exists because 
DomV is semi-open (see the explanations in (P6)). A function f as in (ii) can be 
constructed setting E{t) := the left hand side of the inequality in (3.14). A function 
£ as in (ui) is constructed putting £{r,t) := supj^g^-jy \\V{f,t) — V{fo,t)\\_; this 
sup is proved to be finite using the Lipschitz type inequality (2.14) for V on each set 
C :— B(/o, r) X {t}. One checks that £{r, t) is continuous in (r, t) and non decreasing 
in r. 

Step 2. Proof of (a). Let us pick up any R e (0, R'), and define 

g :[to,T') ^ [0,+oo) , t^ g{t) :^ S{t)+ [ dsu_{t - s)£{R,s) . (3.18) 

J to 

Then Q is continuous and Q{to) = 0; this fact, with the positivity of R, implies the 
existence of T G [to,T') such that Q{t) ^ R for all t G [to,T); the last inequality is 
just the thesis (3.16). 

Step 3. Proof of (b). We apply Proposition 3.7 to the approximate solution 

<^„,(t):=/o for alHG [to,T) . (3.19) 

Due to (i)-(iii), the function S \ [tQ,T) is an integral error estimator for and 
the function £ \ [0,R') x [to,T) is a growth estimator for V from ip^^ (the function 
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p in the general Definition 3.5 of growth estimator is given in this case by pit) — 
constant = i?'). Eq. (3.16) tells us that the general control inequality (3.12) is 
fulfilled by the function Ttit) := constant = R for all t G [to,T). So, Proposition 
3.7 implies the existence of a solution of VP(/o,to)on [to, 7"), and also gives the 
inequality (3.17). □ 



4 Proof of Proposition 3.7. 

Let us make all the assumptions in the statement of the Proposition. We begin 
the proof introducing an appropriate topology for the space of continous functions 
[to,T)^F. 

4.1 Definition. From now on, C{[tQ,T), F) will be viewed as a (Hausdorff, com- 
plete) locally convex space with the topology of uniform convergence on all compact 
suhintervals [to, t] C [to, T). By this, we mean the topology induced by the seminorms 
(II l|T)re[to,T), where 

II ||,:C([to,T),F)^[0,+oo), ^^\\^\\^:= sup (4-1) 

te[to,T] 

(\\ II is the usual norm of F). 

4.2 Remcirk. The uncountable family of seminorms (|| ||^) is topologically equiv- 
alent to the countable subfamiliy (|| ||t-„), where (r„) is any sequence of points of 
[toiT) such that hm„_>+oo t„ = T. Therefore, C{\tQ,T),¥) is a Frechet space. 

To go on we introduce a basic set, whose definition depends on and on the 
function TZ in the control inequality (3.12). 

4.3 Definition. We put 

S) := e C([to,T),F) I VaAm ^ nt) fone [to,T) } . (4.2) 

4.4 Lemma, (i) D is a closed subset of C{[to,T), F), in the topology of Definition 

4.1. 

(a) For all ip eD, one has grip C DomV (and so, J'{ip) is well defined). 

Proof, (i) Suppose ip G C{[tQ,T),F) and ip = lim„^oo'0n, where (ipn) is a sequence 
of elements of D. Then, for all t E [to, T) we have \\ip{t) — ip^p{t)\\ = lim„^oo ll'^n(t) — 

(ii) Let us consider the function p e C([to, T), (0, +oo]) mentioned in the statement 
of Proposition 3.7. Then, t/^ e £) =^ \\ip{t) - </3„j,(t)|| ^ n{t) < p{t) for all t e [to,r) 
^ gYip C 7{ip,^, p) C DomV. □ 
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Prom now on, our attention will be focused on the map 

^^C{[to,T],F) , ^P^Ji^P). (4.3) 

Of course, for (/? e 2), we have the equivalence 

(^solves VP(/o, to) ^ J{^) = ^. (4.4) 

To clarify the sequel, let us recall that a e (0, 1] is the constant appearing in Eq. 
(2.12). 

4.5 Lemma, (i) For each r e \tQ,T) there is a constant A,- e [0, +oo) such that, 
for all e S); 

\\J{m)-Ji^'m\\<K fds forteM . (4.5) 

(ii) For all r e [toiT) and G £), the above equation implies a Lipschitz type 
inequality 

WJi^) - jm\r ^ ^^^^^11^ - v^'ii. (4.6) 

(7 

(which ensures, amongst else, the continuity oj J on D). 

Proof, (i) Let ^,"0' e ®- We consider Eq. (2.23) for J{'ip){t), and subtract from it 
the analogous one for J{'4>'){t). After applying the norm of F and using (2.11), we 
infer 

\\j{m)-jm{m ^ f ci5«_(t-.)i|p(^(.),.)-n^'(.),.)ii_ (4.7) 

for all i e [to-,T). To go on, we fix r e [to, T) and define 

a :={(/,«) eFx[io,r] I ||/-(^„,(s)||^7e(s)}. (4.8) 

This is a closed, bounded subset of F x R, and C T( p) C DomV. Therefore, 
by the Lipschitz property (P6) of V, there is a nonnegative constant L{Ct) = 
such that 

\\V{f,s)-V{f,s)\\_ ^ for (/,s), {f,s) e a ■ (4.9) 

Furthermore, recalling Eq. (2.12) for u_, we see that there is another constant Ur 
such that 

fori'G (0,r-to] . (4.10) 
Inserting Eqs. (4.9) (4.10) into (4.7), we obtain the thesis (4.5) with A^ := UrL^-. 
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(ii) For each t e [to,T], Eq. (4.5) implies 

Jto ~ a 
Taking the sup over t, we obtain the thesis (4.6). □ 

4.6 Lemma (Main consequence of the control inequality for TZ). One has 

J(S)cS). (4.11) 

Proof. Let -0 e 2); then 

J{^^J)-<p^, = [J{ipJ-ipJ + [J{^^;)-J{ipJ] = -E{^J + [J{^^^)-J{<pJ] . (4.12) 

We write this equahty at any t e [to,T), exphcitating J'{ip){t) — J'{(p^p){t); this 
gives 

Jto 

Now, we apply the norm of F to both sides and use Eqs. (3.2) for E{(p^j,){t), (2.11) 
for e^*"*^"^, (3.11) for the growth of V from cp^j,: in this way we obtain 

\\jim)-vaM\^m+ fdsu_it-s)£ms)-^M\\,s) ■ (4.14) 

Jto 

On the other hand, ||-?/'(s)— (/9„p(s)|| ^7l{s) implies £(||-?/'(s) — (/9 <,p(s)||, s) ^ i{7l{s), s); 
inserting this into (4.14), and using the control inequality (3.12) for TZ, we conclude 

\\j{m)-'PaM\^^it)+ [ dsu_{t-s)e{n{s),s)^n{t) , (4.i5) 

Jto 

i.e., Jiij) eD. □ 

The invariance of 2) under JT" is a central result; with the previously shown properties 
of J^, it allows to set up the Peano-Picard iteration and get ultimately a fixed point 
of this map. 

4.7 Definition, (ipk) {k e N) is the sequence in D defined recursively by 

:= Vap , J{Vk-i) {k^l) . (4.16) 
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4.8 Lemma. Let r e [to,T). For all k eN, one has 

y,^^{t)-Mm^^r '^'^rX+'l^^''^ forte[to,r] , (4.17) 
where is the constant of Eq.(4.5) and := max^gj^j -^j £(t). So, 

Ikfe+l - 'PkWr ^ , ©r. Ar r(a) (t - to)" ■ (4-18) 

Proof. Eq.(4.18) is an obvious consequence of (4.17). Wc prove (4.17) by recursion, 
indicating with a subscript k the thesis at a specified order. We have (pi — = 
J{iPap) - = -^(<^ap), whcucc \\^pi{t) " <^o(^)|| ^ ^ (t) ^ S^; this gives (4.17)o. 
Now, we suppose (4.17)^ to hold and infer its analogue of order k + 1. To this 
purpose, we keep in mind Eq. (4.5) and write 

||v'fc+2(t) - <^fc+i(t)|| = \\J{^k+i){t) - J{Vk)m (4.19) 

Jto {t-sy- ^ ''T{ka + l)J,^ {t-sy- 

On the other hand, we have (^) 

i ^'(t-.)---r((fc + i)a + i)^' ' ^'-'"^ 

inserting (4.20) into (4.19) we obtain the thesis (4.17)^+1. □ 

The next (and final) Lemma is a generalization of the inequality (4.18), based on 
the Mittag-Lefiler function (see, e.g., [16]). For any cr > 0, this is the entire 
function defined by 

E.:C-C, ^-^.(^):=Er^^^- (4.21) 
In particular, E„{z) G [1, +oo) for all z G [0, +oo) and Ei{z) = for aU z G C . 



'^To check this, make in the integral the change of variable s = to + x{t — to), with x e [0, 1], 
and then use the general identity 

dxx'^-\l - xf-' = for a,/? > . 

Jo r(Q: + p) 
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4.9 Lemma. For all r e [to, T) and k, k' e N, 

- 'PkWr ^ ^,^^''^M Qrc) , h := min(A;, k') (4.22) 
i [na + i j 

(Qra being defined by (4-18)); this implies that {(fk) is a Cauchy sequence. 
Proof. It suffices to consider the case k' > k (so that h — k). Writing ^p^i — (p^ — 
Yl'j=kifj+'>- ~ fj) "si^g Eq.(4.18) we get 

k'-i 

r{ja+l) 



Vk' - ^kWr <KY1 WZfW^ ■ (4.23) 
j=k 



On the other hand, for each z ^ 0, 



k'—l A +00 ■ +00 



g Tija + 1) ^ § nja +1) ''^^ V{{k + l)a + 1) ^^"^^^ 



z^ 



r(A;(7 + 1) ^ r(£(7 + 1) V{ka + 1)' 

(the last inequality depends on the relation V{a + /3 + 1) ^ V{a + l)r(/3 + 1) for 
^ 0). With z = 9^^, from (4.23) (4.24) we obtain (4.22). This equation, with 
the obvious fact that z^ /T{h(T +1) ^ for /i — > oo and fixed ^ e C, implies 

ll'/'fc' — V'jfcllr — for (/c. A;') — > oo , (4-25) 

for each fixed r e \tQ,T). In conclusion, ((^fe) is a Cauchy sequence. □ 

Proof of Proposition 3.7. {^pk) being a Cauchy sequence, lim^H^oo '^k '■— exists 
in C([io, r), F); (/? belongs to S), because this set is closed. By the continuity of J, 
we have 

J{(p) = lim J{p>k) = lim Lpk+i = p> . (4.26) 

ki-*oo ki-*oo 

Now, recalling (4.4) we get the thesis. □ 

4.10 Remark. The Mittag-Leffler function E^- on [0, +oo) is strictly related to a 
linear integral equation. More precisely, given cr > let us consider the following 
problem: find G e C([0, +oo), R) such that 
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This has a unique solution 



G{t) := E^it") for t e [0, +00) . (4.28) 

One checks directly that the above G solves (4.27) (^); uniqueness of the solution 
follows from the linearity of the problem and from Lemma 2.6. 
Integral equations like (4.27) are related to the so-called "fractional differential equa- 
tions" (see e.g. [4], also mentioning E^^). 

5 Applications of Proposition 3.7 to systems with 
quadratic nonlinearity: local and global results. 

The setting. Throughout this section we consider a set {F_^.,F,F_, A,u,u_,'y,^) 
with the following features. 

F+,F,F_ are Banach spaces, A is an operator and u,u_ are semigroup estimators 
fulfilling conditions (P1)-(P5). Furthermore: 
(Ql) y is a bilinear map such that 

yrFxF^F., {f,9)^nf,9); (5-1) 

we assume continuity of which is equivalent to the existence of a constant K e 
(0, +00) such that 

\\m9)\\-<K\\f\\\\g\\ (5.2) 

for all e F (6). 
(Q2) We have 

eeC°'^([0,+oo),F_) , (5.3) 
(recall that C°'^ stands for the locally Lipschitz maps). 
Having made the above assumptions, let us fix some notations. 

5.1 Definition. From now on: 

(i) U e (^([0, -l-oo), [0, +00)) is a nondecreasing function such that 

[ dsu_{s) ^U{t) fort e[0, +00) , U{0)=0 (5.4) 
Jo 

^To this purpose one inserts into (4.27) the series expansion coming from (4.21), and then uses 
the identity in the previous footnote. 

^Of course, in the trivial case 3* = we could fulfill (5.2) with K = as well. In the sequel we 
always assume ii' > 0, to avoid tedious specifications in many subsequent statements and formulas. 
In any case, such statements and formulas could be extended to = by elementary limiting 
procedures. 
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(e.g., U{t) := j^dsu_{s); in any case, (5.4) and the positivity of u_ imply U(t) > 

for allt > O;. We putU{+oo) := limt_+oo G (0,+oo]. 

(ii) H_ e C([0, +00), [0, +00)) is any nondecreasing function such that 

\m\\-<^-it) fort e[0,+oo) (5.5) 

(e.g., E_{t) := sup^gjo^t] ||C(s)|Uy'. We put E:_(+oo) := lim^^^+oo 2_(t) G [0, +00]. 
(Hi) V is the (affine) quadratic map induced by T and ^ in the following way: 

V:Fx[0, +00) ^ F_ , if, t) ^ r{f, t) := /) + ^ ■ (5-6) 

Properties of 'P. Let us analyse this map, so as to matcli tlie scliemes of tlie 
previous sections. First of all we note that DomV is semiopen in F x R, in the 
sense defined in (P6). Furtheremore, we have the following. 

5.2 Proposition. For all {f,t) and {f',t') e Fx [0,+oo), 

\\V{f, t) - V{f\ t')\\_ ^ 2K\\f\\ 11/ - /'II + K\\f - ff + \m - i{t')\\_ . (5.7) 
Proof. For the sake of brevity, we define h :— f — f . Then 

v{f,t) = yu' + Kf' + h) + i{t) 

= ?(/', /') + ?(/', h) + f) + nK h) + m ■ 

subtracting V{f' ,t'), we get 

vu. t) - v{f', t') = nf, h) + f) + k) + m - m . 

We apply || ||_ to both sides, taking into account Eq. (5.2); this gives \\P{f,t) — 
V{f\t')\\_ < 2ii'||/'|| \\h\\ + K\\hf + Wiit) - yielding the thesis (5.7). □ 

The previous proposition has two straightforward consequences. 

5.3 Corollary. For each hounded set C of F x [0, +00), there are two constants 
L = L{C), M = M{C) such that 

\\V{f, t) - Vif, t')\\_ ^L\\f-f'\\+M\t- t'\ for (/, t), (/', t')eC; (5.8) 

so, V fulfills condition (P6). Let us denote with B and I the projections of C on F 
and [0, +00), respectively. Then we can take L :=4i^||i3||, t(;/iere ||i3|| :=supjgg||/||; 
furthermore, we can take for M any constant such that \\iit) — ^{t')\\_ ^ M\t — t'\ 
for t,t' e I. 

Proof. In Eq. (5.7), we substitute the relations ||/'|| ^ ||i3||, ||/ - /'|p ^ (||/|| + 
11/11)11/ - /II ^ 211^1111/ - /'II, and the inequality defining M. □ 
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5.4 Corollary. Let us consider any function (p e C{[0,T), F) . Then 
\\V{f,t)-Vm),t)\\_ ^ £{\\f-(l>{t)\\,t) for aUfeF,te [0,T) , (5.9) 

£ : [0, +00) X [0, T) ^ [0, +00) , (r,t)^£(r,t):^2K\\(l)(t)\\r + Kr'^ . (5.10) 

The function i is a growth estimator for V from (f), in the sense of Definition 3.5 
(with a radius of the tube p{t) :— +00 for all t). 

Proof. Use Eq. (5.7) with (/',t') := {(j){t),t). □ 

Cauchy and Volterra problems; approximate solutions. These problems will 
always be considered taking to := as the initial time; we will write 

VP(/o) := VP(/o, 0) , CP(/o) := CP(/o, 0) for each /o . (5.11) 

For the above problems, we have the following results. 

(a) If /o e F+ and F_ is reflexive, VP(/o) is equivalent to CP(/o) (see Proposition 
2.4). 

(b) For any /o G F, uniqueness and local existence are granted for VP(/o) : see 
Propositions 2.5 and 3.10. 

(c) We can apply to VP(/o) Proposition 3.7 on approximate solutions, choosing ar- 
bitrarily the approximate solution (f^^; as an error estimator, we can use at will the 
function £ in Corollary 5.4 (or any upper bound for it). This yields the following 
statement. 

5.5 Proposition. Let fo^F and T e (0,-|-oo]. Let us consider for VP(/o) an 

approximate solution (p^p e C([0,T), F), and suppose there are functions £,T>,TZ & 

C([0,T), [0,+oo)) such that (i)-(ui) hold: 
(i) ipap hcLS the integral error estimate 

\\E{^^M\\<^it) forte[0,T) ; (5.12) 

(a) one has 

||V9„,(t)Kl)(t) forte[0,T); (5.13) 

(Hi) with K as in (5.2) andU as in (5.4), solves the control inequality 

S{t) + K [ dsu_{t- s){2V{s)n{s) +n\s)) ^n{t) forte[0,T). (5.14) 
Jo 

Then, (a) and (b) hold: 

(a) VP(/o) has a solution (p : [0, T) F; 

(b) one has 

\\ip{t)-ip^M\^nt) forte[0,T) . (5.15) 
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Proof. We refer to the control inequality (3.12) in Proposition 3.7 (with to — 0). 
Due to Corollary 5.4, the growth of V from ip^^ has the quadratic estimator i{r, t) :— 
2K\\(p^p{t)\\ r + Kr'^; binding ||</7ap(^)|| via (5.13) we get another estimator, that we 
call again i, of the form 

i{r, t) = 2KV{t)r + Kr^ . (5.16) 

With this choice of £, the control inequality (3.12) takes the form (5.14) and (a) (b) 
follow from Proposition 3.7. □ 

How to handle the control inequality (5.14). Rephrasing Remark 3.9 in the 
present case, we repeat that TZ is the only unkown in (5.14). In fact, the functions 
S,V appearing therein can be determined when ip^^ is given, and u_, K can be 
obtained from A, ^ (as an example, the computation of u_ and K for the NS 
equations will be presented in Sections 6-7) . Two basic strategies to find a function 
IZ solving (5.14) on an interval [0, T), if it exists, can be introduced: 

(a) the analytical strategy: one makes some ansatz for 7?., substitutes it into (5.14) 
and checks whether the inequality is fulfilled; 

(b) the numerical strategy. 

Let us mention that a numerical approach was presented in [12], for the simpler 
control inequality considered therein; in that case it was possible to transform the 
control equality (with ^ replaced by =) into an equivalent Cauchy problem for TZ, 
and then solve it by a standard package for ODEs. 

The approach of [12] can not be used for (5.14) due to the singularity of u_{t) for 
f — > 0"^; a different numerical attack could be used, but this is not so simple and 

its features suggest to treat it extensively elsewhere. For the above reasons, in the 
present work we only give an introductory sketch of the numerical strategy: see 
Appendix B. In the rest of the paper, starting from the next paragraph, we will use 
the analytical strategy (a). 

Some special results on VP(/o) . All these results will be derived solving the 
control inequality (5.14) by analytic means, in special cases. 

5.6 Proposition. Let /q E F and T E (0,+oo]. Let us consider for VP(/o) an 
approximate solution (f^p E C{[0,T), F), and suppose there are functions £,V E 
C([0,T], [0, +oo)) such that (i)-(iii) hold: 

(i) £ is nondecreasing and binds the integral error as in (5.12); 

(a) V is nondecreasing and binds ip^^ as in (5.13); 
(Hi) with K as in (5.2) andU as in (5.4) C) ; 

2^KU{T)£{T) + 2KU{T)V(T) ^ 1 . (5.17) 
''Of course, in the case T = +oo (5.17) implies U{+(X)) < +oo 
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Then VP(/o) has a solution ip : [0,T) ^ F and, for all t e [0, T), 

Mt)-^^M\<nt) , (5-18) 

f l-2KUit)Vit) - ^{l-2KU{t)V{t)Y-AKU{t)S{t) 

U(0) ift^Q; 

the above prescription gives a well defined, nondecreasing function TZ e C{[0,T), 
[0,+oo)). 

Proof. We refer to Proposition 5.5, and try to fulfill the control inequahty (5.14) 
with a nondecreasing U e C([0, T), [0, +oo)). Noting that 7^(s) ^ 7^(^), V{s) ^ 
'D{t) for s e [0,t], we have 

£{t) + K f dsu_(t- s){2V{s)n{s) +n'^{s)) (5.19) 
Jo 

^ £{t) + K{2V{t)TZ{t) +n'^{t)) [ dsu_{t-s) 

Jo 

^ £{t) + K{2V{t)n{t) + 7^2(^)) u{t) 

the last inequality follows from j^dsu_{t — s) = dsu_{s) ^ W(t). Due to (5.19), 
(5.14) holds ii £{t) + K{2U{t)V{t)n{t) +U{t)n'^{t)) ^ 7^(^), i.e., 

KU{t)n{tf - (1 - 2KU{t)v{t))n{t) + £:(i) ^ o . (5.20) 

This inequality is fulfilled as an equality if we define 7?. as in (5.18), provided that 
2y/KU{t)£{t) + 2KU{t)V{t) ^ 1; this happens for each t e [0,7) due to the as- 
sumption (5.17) (^). 

The function TZ on [0,T) defined by (5.18) is continuous and nonnegativc; to con- 
clude the proof, we must check it to be nondecreasing. To this purpose, we note 
that 

n{t) = T{KU{t),V{t),£{t)), (5.21) 



r 1-2/1(5- V(l- 2^1(5)^ -4/X6 

[ e if /X = 0. 

The above function T has domain 

DomT := {(/x, 6, e)\ ii,6,e^ 0, 2^6 + 2^ ^ 1 } , (5.22) 



^Obviously enough, we take for TZ{t) the definition (5.18) since this gives the smallest nonneg- 
ative solution of Eq. (5.20). 
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containing all triples {KU{t),V{t),£{t)) due to (5.17); one checks by elementary 
means (e.g., computing derivatives) that T is a nondecreasing function of each one 
of the variables /i, S, e, when the other two are fixed. □ 

5.7 Remark. The inequality (5.17) is certainly fulfilled if {S{T),T>[T)) or T are 
sufficiently small (recall that U{T) vanishes for T — > 0^). 

An example: the zero approximate solution. For simphcity, we suppose 

u{t) ^ 1 for all t e [0, +oo) . (5.23) 

Let /o G F, T G (0, +oo]; we apply Proposition 5.6, choosing for VP(/o)the trivial 
approximate solution 

ipUt):=0 forie[0,T). (5.24) 

5.8 Lemma, (p^p '■— has the integral error 

E{^J{t) = -e'^fo - f dse-^'-^^^i{s) ; (5.25) 

with H_ as in Definition 5.5, ||£'((/9„p)(t) || has the estimate 

\\E{ipJ{t)\\ ^ J{t) , ^{t) \\h\\ + ^_{t)U{t) . (5.26) 

Proof. Eq. (5.25) follows from the general definition (3.1) of integral error, and 
from the fact that V{(pap{s), s) = ^{s). 

Having Eq. (5.25), we derive (5.26) in the following way. First of all, 

\\E{^jm<u{t)\\M\+ f ds u_{t - s)U{s)\\_ ; (5.27) 
but u{t) ^ 1, U{s)\\_ ^ E_{t) for s G [0,t], so 

\\E{vJ{t)\\^\\fo\\+^-{t) fdsu_{t-s)^\\fo\\+E_{t)U{t). □ (5.28) 

Jo 

From the previous Lemma and Proposition 5.6, we infer the following result. 

5.9 Proposition. With u, 3^ as in (5.23) (5.26), assume 

AKU{T)J{T) ^ 1 . (5.29) 
Then VP(/o) has a solution (p : [0, T) ^ F and, for all t in this interval, 

\\Lp{t)\\ ^ "Jit) X{AKU{t)-J{t)) . (5.30) 
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Here A" e C([0, 1], [1, 2]) is the increasing function defined by 

X{z) := I ' "(./2)" ' /^^^^(O'l]' (5.31) 
[ 1 for z — . 

Proof. According to the previous Lemma, we can apply Proposition 5.6 with 99 = 
0, S = 9^; obviously, we have for ||y9„p(t)|| the estimator V^t) := 0. Eqs. (5.17), 
(5.18) yield the present relations (5.29), (5.30), (5.31); in particular, the function TZ 
in (5.18) is given by 7^(^) = - iKU{t)3^{t)) /2KU{t) = 3^(t) X{AKU{t)'J{t)). 

□ 

Of course, the basic inequality (5.29) is fulfilled if (/o,S_(T')) or T are sufficiently 
small. 

Further results (global in time) for VP(/o) , under special Eissumptions. 

We keep the assumptions (P1)-(P5) of Section 2 and (Q1)(Q2) at the beginning of 
this section, and put more specific requirements on the semigroup estimators u, u_. 
More precisely, we add to (P4) (P5) the following conditions, involving two constants 

S e [0, +00) , AT e (0, +00) . (5.32) 
(P4') The semigroup estimator u has the form 

u{t) = e"-^* for t e [0, +00) . (5.33) 
(P5') The semigroup estimator u_ has the form 

u_{t) = fx_{t) e-^* for t e (0, +00) , (5.34) 

//_eC((0,+oo),(0,+oo)), /x_(t) = 0(^) fort^O+ (a e [0, 1)) , (5.35) 

f ds ^x_{t - s) e-^' for t e [0, +00) . (5.36) 

Jo 

5.10 Proposition. Given /q G F, let us consider for W{fo) an approximate solu- 
tion ip^p G C([0, +00), F). Suppose there are constants E,D E [0, +00) such that: 
(i) admits the integral error estimate 

||E((^J(t)|| ^ ^e-^* forte [0,+oo) ; (5.37) 

(a) for all t as above, 

y^M\<De-^' ; (5.38) 
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(Hi) with N as in (5.36) and K as in (5.2), one has 

2VKNE + 2KND ^ 1 . (5.39) 
Then VP(/o) has a global solution (p : [0, +oo) F and, for all t e [0, +oo); 

Mt) - V.M\ < Re-""' : (5.40) 



_ 1 - 2KND - ./(^I^^WNDY^^IKNE 
'~ 2KN ■ 

Proof. Step 1: the control inequality. We use again Proposition 5.5 and the 
control inequality (5.14). With the notations of the cited proposition, wc have 
S{t) := Ee-^\ V{t) := De"-^* and u_ has the expression (5.34). So, (5.14) takes 
the form 

Ee-^* + K f ds ix_{t - s)e-''^'-'\2De~'''n{s) +Te{s)) ^n{t) , (5.41) 
Jo 

that we regard as an inequality for an unknown nonnegative function 71. 
Step 2: searching for a global solution TZ of (5.41). We try to fulfill (5.41) with 

7^(^) Re-^^ for all t e [0, +oo), (5.42) 

with an unknown constant. Then, the left hand side of (5.41) is 



E + K{2DR + R^) f dsfi_{t 
Jo 



s)e 



(5.43) 



^ e"^* [E + K{2DR + R'^)N) , 



where the inequality depends on (5.36). The last expression is bounded by TZ{t) if 
R fulfills the inequality E + K{2DR + R^)N ^ R, i.e., 

KNR^ - (1 - 2KND)R + ^ . (5.44) 

This condition is fulfilled as an equality if we define R as in (5.40); this R is well 
defined and nonnegative due to the assumption (5.39) (^). Due to the above con- 
siderations, the thesis is proved. □ 

Applications to cases with exponentially decaying forcing. 

From here to the end of the section, we add to (P1)-(P5')(Q1)(Q2) the following 
condition: 



^and, in fact, is the smallest nonnegative solution of (5.44). 
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(Q3) There is a constant J e [0, +00) such that (with B as in (P4')(P5')) 

U{t)\\. ^ for all t e [0, +00) . (5.45) 

Two cases where (Q3) holds are: (i) the trivial case ^{t) — for all t; (ii) situations 
where the external forcing is switched off in the future. 

Hereafter we present two applications of Proposition 5.10, corresponding to different 
choices for Both of them give global existence for the exact solution (/? of 

VP(/o) when the datum /o is sufficiently small, with suitable estimates of the form 
(5.40). 

Example: the zero approximate solution. Let /o e F; we reconsider, from the 
viewpoint of Proposition 5.10, the VP(/o) approximate solution 

^,^{t) := for t e [0, +00) . (5.46) 

5.11 Lemma, (f^p :— has the integral error estimator 

\\E{ipJ{t)\\ ^ Fe-^' , F := ||/o|| + NJ . (5.47) 

Proof. The integral error E{ip^p) was already computed, see Eq. (5.25). Prom this 
equation and the assumptions (5.33) (5.34) on u, u_ we infer 

\\E{ipJ{t)\\^e-^'m+ fdse-^^'-^^i,_{t-sM{s)\\_ ; (5.48) 
inserting here the bound (5.45) for ||(^(s)||_, and using Eq. (5.36) for /x_ we get 

\\E{^J{t)\\^e~'''\\fo\\ + Je-^' f ds i^_{t-s)e-''' ^e-'^'WhW + Je-^'N . (5.49) 

Jo 

□ 

Prom the previous Lemma and Proposition 5.10, we infer the following. 

5.12 Proposition. With F as in (5.Ji.'7), let 

AKNF ^ 1 . (5.50) 
Then VP(/o) has a global solution ip : [0, +00) F and, for all t e [0, +00), 

y(t)\\ ^ FX{AKNF) e'^* , (5.51) 

with X as in (5.31). 
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Proof. We apply Proposition 5.10 with (p^^ :— 0. The constants of the cited 
proposition are E — F and D — Q (the first equaUty follows from the previous 
lemma, the second one is obvious). Eqs. (5.39) (5.40) yield the present relations 
(5.50) (5.51) (5.31); in particular, the constant in (5.40) is given by = (1 — 

VI -AKNF)/{2KN) = FX{AKNF). □ 

Example: the ".4, -flow" approximate solution. Given /o € F, we consider 
for VP(/o) the approximate solution 

ip^^it) := e*^/o + /* ds e(*-^)^^(s) for t e [0, +oo) (5.52) 

(i.e., we use the flow of the linear equation / = Af + ^{t}). 

5.13 Lemma. Let F be as in (5.47). The above ip^p has the integral error estima- 
tor 

\\E{^J{t)\\^KNF^e-^' (5.53) 
and fulfills for all t e [0, +oo) the norm estimate 

\\ip^M\^Fe-^' ■ (5-54) 

Proof. We first derive Eq. (5.54). To this purpose, we note that the definition 
(5.52) of (pap, the assumptions (5.33) (5.34) on u,u_ and Eqs. (5.45) for ||^(-)||_, 
(5.36) for /x_ imply 

yapim ^ e-^*||/o|| + tds e-^(*-)/._(t-.)||e(.)||_ (5.55) 

^0 

^ e-^*ll/o|| + Je'""' f ds fi_{t - s)e-^^ ^ e'^'^ll/oll + Je-'^'N ; 
Jo 

by comparison with the definition (5.47) of F, we get the thesis (5.54). 

Let us pass to the proof of (5.53). To this purpose we note that the definition (3.1) 

of integral error gives, in the present case, 

E{ipJ{t) = - f dse^'-'^-^'y{^^,{s),pM) ■ (5.56) 
Prom here and from Eqs. (5.34) for m_, (5.2) for (P we infer 

\\E{Vap)m ^ fds i,_{t- s)e-^^'-%nvM.<Pap{s))\\- (5.57) 

fds ^_(t-5)e-^(*-^)||^.,(s)||2 . 
Jo 
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In the last inequality, we insert the bound (5.54) and then recall (5.36). This gives 

\\E(^J (t) II ^ KF^e-^' f ds e-^'ii_{t - s) ^ RNF^e'^' . □ (5.58) 

Jo 

Let us return to Proposition 5.10; with the previous Lemma, this implies the fol- 
lowing result. 

5.14 Proposition. Let us keep the definition (5.47) for F, and the assumption 
(5.50) AKNF ^1. The global solution cp : [0, +oo) F o/ VP(/o) is such that, for 
all t e [0,+oo), 

Wifiit) - e*-^/o|| ^ KNF^ XiAKNF) e"^* ; (5.59) 
here X e C([0, 1], [1,4]) is the increasing function defined by 

[ 1 forz^O. 

Proof. According to the previous Lemma, we can apply Proposition 5.10 with 
E = KNF^ and D^F. Eqs. (5.39) (5.40) yield the present relations (5.50) (5.59) 
(5.60); in particular, the constant in Eq. (5.40) is 

R={1- 2KNF - Vl - AKNF)/{2KN) = KNF^ X{AKNF) . □ 

5.15 Remark. Most of the results presented in this section could be extended to 
the case P(/, t) = 7{f, f)-\-^(t), where J" : F™ F_ is continuous and m-linear 
for some integer m ^ 3. In this case, the growth of V from any approximate solution 
admits an estimator £{r, t) more general than (5.10), which is polynomial of degree m 
in r. One could extend the analysis as well to the case V{f, t) = ?"(/, /, t) +C,{t), 
involving a time dependent multilinear map T : F™ x [0, +oo) — » F_, (/i, fm, t) i— > 
y(/i, fm,t). These generalizations are not written only to save space. 

6 The Navier-Stokes (NS) equations on a torus. 

Prom here to the end of the paper, we work in any space dimension 

d^2 . (6.1) 

Preliminciries: distributions on T*^, Fourier series and Sobolev spaces. 

Throughout this section, we use r, s as indices running from 1 to d and employ for 
them the Einstein summation convention on repeated, upper and lower indices; 5rs 
or S^^ is the Kronecker symbol. 
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Elements a,b,.. of R*^ or will be written with upper or lower indices, according 
to convenience: (a*) or (og), (b^) or {bg). In any case, a»b is the sum of product of 

the components of a and b, that will be written in different ways to accomplish with 
the Einstein convention. Two examples, corresponding to different positions for the 
indices of a, are 

a»b — Qgb^ , a»b — Srsci^b^ ■ (6.2) 
Let us consider the (i-dimensional torus 

T'^ := T X ... X T {d times) , T := R/(27rZ) . (6.3) 

A point of T*^ will be generally written x — {x^)r=i,...d- We also consider the "dual" 
lattice Z'' of elements k — {kr)r=i,...,d and the Fourier basis {ekjk^^d, made of the 
functions 

e.:T'^-C, e.(a;):=^e^^- (6.4) 

{k»x := krX^ makes sense as an element of T'^). We introduce the space of periodic 
distributions D'iT"^, C) = D'^, which is the dual of C°°(T'^, C) = (equipping 
the latter with the topology of uniform convergence of all derivatives); we write 

(f , /) for the action of a distribution v G D'q on a test function / G C^. The 
weak topology on Dq is the one induced by the seminorms pf {f & Cq), where 
Pfiv) ■■=\{vJ)\. 

Each V e D'q has a unique (weakly convergent) series expansion 



VkCk , (6.5) 



kezd- 

with coefficients Vk & C for all k, given by 

Vk = {v,e_k) ■ (6.6) 

The "Fourier series transformation" v ^ {vk) is one-to-one between Dq and the 
space of sequences s'(Z'^, C) = s'^., where 

s'q :— {c — (cfe)fe£zd I Cfe e C, |cfc| = 0(|/c|^) as A; — > oo, for some p G R} . (6.7) 

In the sequel, we often use the mean of a distribution v G D'^, which is 

(in the first passage, {v, 1) means the action of v on the test function 1; the second 
relation follows from (6.6) with k = 0, noting that cq — 1/(277)"^/^. Of course, 
(v, 1) = v{x)dx if V is an ordinary, integrable function). 
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The complex conjugate of a distribution v e D'(^ is the unique distribution v such 
that {v, f) = {v, /) for each / e Cg"; one has = J^kez^ ^k^-k- 
From now on we will be mainly interested in the space of real distributions D'(T'^, R) e 
D', defined as follows: 

D' :^{veD'c\v^v}^{veD'c\v^^ v^k for all k e Z'^} ; (6.9) 

we note that v & D' implies (v) £ R. The weak topology on D' is the one inherited 
from Dq. 

Let us write ds {s — l,...,d) for the distributional derivative with respect to the 
coordinate x^; from these derivatives, we define the distributional Laplacian A := 
6^^drds : D'fj D'fj. Of course dgCk — ikgCk, Ae^ = —\k\'^ek for each k. For any 
V e Dq, this implies 

dgv^i^ ksVkek , Av ^ - ^ \k\hkek ; (6.10) 

{1-Arv= J2{l + \k\'rvkek (6.11) 
fceT<* 

for m e 0,1,2,3,.... For any m e R, we will regard (6.11) as the definition of 
(1 — A)"* as a linear operator from D'(~. into itself. Comparing the previous Fourier 
expansions with (6.8), we find 

{dsv) = , {Av) = . (6.12) 

All the above differential operators leave invariant the space of real distributions, 
more interesting for us; in the sequel we will fix the attention on the maps ds, A, 
(1 - A)™ : D' D'. 

Let us consider the real Hilbert space L^(T'^, R, dx) = L^, i.e., 
L'^ :^ {v: T'^ ^ R \ [ v'^{x)dx < +00} ^ {v e D' \ l^^l^ < 5 (6-13) 

this has the inner product and the associated norm 



{v\w)l2 / v{x)w(x)dx — 2_] VkWk , (6-14) 



\ v'^{x)dx ^ \vk\'^ . (6.15) 
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To go on, we introduce the Sobolev spaces if"(T'^, R) = if". For each n e R, 
ii" := {veD'l (l-A)"/2^ eL^}^{veD'\ J2 (l + I^D^I^fel' < +00 } ; (6.16) 

this is also a real Hilbert space with the inner product 

{v\w)^ := ((1 - Ar/\ I (1 - Ar/'w)L2 = J] (1 + \k\'r-kw, (6.17) 

and the corresponding norm 



||^;|U := 11(1 - ^lU. = J E (1 + ■ (6-18) 

One proves that 

n^n' ^ H^^H""' , \\ ^ || ||„ . (6.19) 

In particular, is the space and contains if" for each n ^ 0; for any real n, A 
is a continuous map of into H"~^. Finally, let us recall that 

//" ^ D' for each n e R ; (6.20) 

H^'^C'^ ifgeN, n e (g + +oo) . (6.21) 

In the above if" carries its Hilbertian topology, and D' the weak topology; C stands 
for the space C"'(T'^, R), with the topology of uniform convergence of all derivatives 
up to order q. 

Obviously enough, we could define as well the complex Hilbert spaces and ifg; 
however, these are never needed in the sequel. 

Spaces of vector valued functions on T*^. To deal with the NS equations, we 
need vector extensions of all the above spaces and mappings. Let us stipulate the 
following: if V{T'^, R) = is any vector space of real functions or distributions on 
T"^, then 

Y:=V'^={v = {v\...,v'^)\v' eV for all r}. (6.22) 

This notation allows to define the spaces D', L^, H". Any v = (v^) G D' will be 
referred to as a vector field on T''. We note that v has a unique Fourier series 
expansion (6.5) with coefficients 

Vk = {vl)r=,,...,d e C^ vl := {v\ e.k) ; (6.23) 

again, the reality of v ensures Uk = v^k- We define componentwisely the mean 
{v) e R'^ of any v e D' (see Eq. (6.8)), the derivative operators dg : li' —>■ D', their 
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iterates and, consequently, the Laplacian A. The prescription (6.11) gives a map 
(1— A)"* : B' — > D' for all real m. Whenever V is made of ordinary functions, a d-uple 
V G V can be identified with a function v : ^ H'^, x i— v{x) = {v^'{x))r=i,...^d- 

is a real Hilbert space. Its inner product is as in (6.14), with v{x)w{x) and VkWk 
replaced by 

v{x)»w{x) — SrsV'^{x) W^{x) , Vk»Wk — SrsVlwl ; (6-24) 

the corresponding norm is as in (6.15), replacing v^{x) with |f(a;)p = Ylt=i'^^ i-'^)^ 
and intending \vk\^ = J2t=i I'^fcP- 

For any real n, the Sobolev space H" is made of all d-uples v with components 

G H"", an equivalent definition can be given via Eq.(6.16), replacing therein 
with L^. is a real Hilbert space with the inner product 

{v\w)n := ((1 - A)"/2^ 1(1 - A)"/2^)i. = E (1 + 1^1')"^-^^ ■ (6-25) 

The corresponding norm || ||„ is given, verbatim, by Eq. (6.18); Eq. (6.19) holds as 
well for H", H'* and their norms. Let us consider the Laplacian operator A : D ^ D ; 
for any real n 

AH" C H"-2 , (6.26) 

and A is continuous with respect to the norms || ||„, || \\n-2- The embeddings (6.19) 
(6.20) (6.21) have obvious vector analogues. 

Zero mean vector fields. The space of these vector fields is 

:= G D' I (v) = 0} (6.27) 

(of course, (v) = is equivalent to the vanishing of the Fourier coefficient Vq). 

Divergence free vector fields. Let us consider the divergence operator (hnear, 
weakly continuous) 

div :W^D', v^diw.^drv'' ; (6.28) 

we put 

D'j, := {v eW \ divv = 0} (6.29) 

and refer to this as to the space of divergence free (or solenoidal) vector fields. The 
description of these objects in terms of Fourier transform is obvious, namely: 

divv — i {k»Vk)ek for all v — VkCk G D' ; (6.30) 
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= {v eB' \ k.Vk = for all k e Z'^ } (6.31) 
^ {v \ Vk e-< k for all A; e Z'^ } , 

where -< k y is the subspace of spanned by k, and is the orthogonal comple- 
ment with respect to the inner product {b, c) G x 6» c. 
In the sequel, we will consider as well the subspace 

D',„:=D',nD^ (6.32) 

Gradient vector fields. Let us consider the gradient operator (again linear, weakly 
continuous) 

d:D'^B>' ; p ^ dp := idsP)s=i,...,d ; (6.33) 
if P = ^keziPk^k, then 

dp = i ^ kpkCk . (6.34) 
kezd- 

The image 

D; ^{dp\pe D'} (6.35) 

is a linear subspace of D', hereafter referred to as the space of gradient vector fields; 
for any vector field w, comparison with (6.34) gives 

B>'^^ {w eB>' \ Wk e< k y for all k e Z'^} . (6.36) 

Of course, if w is in this subspace, the distribution p such that w — dp is defined up 
to an additive constant. We put 

d~^w :— unique p & D' such that w — dp and po = . (6.37) 

This gives a linear map 

d-^ : D; ^ D' . (6.38) 

The Leray projection. Using the Fourier representations (6.31) (6.36), one easily 
proves the following facts. 

(i) One has 

©' = D'j. © ©; (6.39) 

in algebraic sense, i.e., any v e D' has a unique decomposition as the sum of a 
divergence free and a gradient vector field. 

(ii) The projection 

£ : D' ^ D'^ , v^2v (6.40) 
corresponding to the decomposition (6.39) is given by 
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£^v^J2 i^kVk)ek for all t; = ^ VkCk e D , (6.41) 

£fc := orthogonal projection of onto -< k y-^ ; 
more explicitly, for all c e C^, 

£oC = c, £fcC = c-^^^ for A; e Z'^,^ ^ . (6.42) 

As usually, we refer to £ as to the Leray projection; this operator is weakly contin- 
uous. From the Fourier representations of £, of the mean and of the derivatives one 
easily infers, for all v e D', 

(£^;) = (v) , £{dsv) = d,{2v) , £(A^;) = A{£v) . (6.43) 

A Sobolev framework for the previous decomposition. For n e R, let us 

define 

:= H" n D'j, = e M" I divt; = 0} ; (6.44) 

H;! IT n = {w; G H" | = dp, p e D'} ; (6.45) 

ai/" := {ap I p e ii"" } . (6.46) 
Then the following holds for each n: 

(i) Hp is a closed subspace of the Hilbert space (H", ( | )„) (because div is continuous 
between this Hilbert space and D' with the weak topology). 

(ii) One has 

= 9//"+^ (6.47) 

and H" is also a closed subspace of H". The map of Eq. (6.37) is continuous 
between and //"+^ 

(iii) Denoting with the orthogonal complement in (H", ( | )„), we have 

and £ \ H" is the orthogonal projection of H" onto H^; so, as usual for Hilbertian 
orthogonal projections, 

ll-C^^lln ^ Iklln for all V eW . (6.49) 

Other Sobolev spaces of vector fields. For n e R, we put 

:= W n D; := {/ G H" | (/) = 0} ; (6.50) 

mi, := mi n m^ -.= {/ e | div/ = o , (/) = 0} (6.51) 
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Then, is a closed subspace of H" (by the continuyity of ( ) : H" — > C) ; the same 

holds for H^g, since this is the intersection of two closed subspaces. 

The space (6.51) plays an important role in the sequel; we will often use the Fourier 

representation 

= {/ e » I E (1 + < +~ ' (6-52) 

fc./fc = for all k E Z'^Jq = 0} . 

Some inclusions. We note that the relations AH" C HI"~^ and div(Av) = A(divi'), 
(Av) = for all v e D imply the following, for each real n: 

AH" c H^-^ , AH^ c H^-^ . (6.53) 

A digression: estimates on certain series. Let us define 

- \ {0} ; (6.54) 
throughout the section, n is a real number such that 

n > ^ . (6.55) 

The series considered hereafter are used shortly afterwards to derive quantitative 
estimates on the fundamental bilinear map appearing in the NS equations: by this 
we mean the map sending two vector fields v,w on T'' into the vector field v» dw 
(see the next paragraph). The estimates we give are also useful for the numerical 
computation of those series. In both Lemmas hereafter, 

Z'^ := Z'^ or . (6.56) 

6.1 Lemma. One has 

For any real "cutoff" X ^ 2y/d, one has 

S„(A) < S„ ^ S„(A) + SSnW (6.58) 

where 

5BM := 0^±dy^ _ .ggQX 
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Proof. See Appendix C. □ 
6.2 Lemma. For k e Z'^, define 

then, (i) (ii) hold. 

(i) One has K,n{k) < +00 for all k e Z"^; furthermore, with E„ as in (6.57), 

Kn{k) —>■ S„ for /c — > 00 . (6.62) 

Thus, 

sup )Cn{k) < +00 . (6.63) 

(ii) Let us choose any " cutoff function" A : Z^ ^ [2-\/d, +00) and define 

( i + \k\' ifm<\k\ , 

A : ^ (0, +00) , k^X{k):={ 1 + \k\' if A(k) > \k\ ^^"^^^ 

i l + (A(fc)-|A;|)2 ^^^^^^^1^1 ■ 

r/ien, for all k e Z"^, 

Xnik) < lCn{k) ^ DC„(A;) + 5DC„(A;) , (6.65) 

X(,y-(l±mi V l'^-''!' (666) 

53C„(*) := (l±irmrl . (6.67) 

^ 2d-in^/2r{d/2){2n-d) {A{k)-Vdy^-d ^ ' 

Finally, suppose the cutoff A has the property 

a\k\ ^ K{k) ^ I3\k\ for all keZ'^ such that \k\ ^ x (6-68) 
(X ^ 0, 1< a ^ /?) ; 



3<;„(A;) ^ S„ , 53C,(A;) = C'(^^^) ^0 /or A: ^ 00 . (6.69) 
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Proof. See Appendix C. □ 

The fundamental bilinear map. Let v,w E H". For each r,s E {1,2, d}, 
E and drW'^ E H^~^ are ordinary real functions: note that v'^ E by the 
embedding (6.21), and dj-w^ E since n — l> d/2 — 1^0. These functions can 
be multiplied pointwisely, and this allows to define 

V, dw := {v. dw')^^^^^^^^^ v. dw' := v'drw' : ^ C . (6.70) 

6.3 Proposition, (i) Consider v,w E H". The vector field v»dw has Fourier 

coefficients 

{v'dw)k^ j^;^ J2[i'"h'{k-h)]wk-h ■ (6.71) 

Furthermore, v»dw E H""^. 
(a) The bilinear map 

WxW^W-^, {v,w)^v.dw (6.72) 
admits an estimate (indicating continuity) 

\\v»dw\\n-l ^ Kn \\v\\n\\w\\n (6.73) 

for all v,w as above, with a suitable constant = Knd G (0, +00). For the latter 
one can take any constant such that 



sup }Cn{k) ^ Kn , (6.74) 

i2^r ,^.(i + iwi+i^-/^iT ^ ^ ^ 

(as in (6.61), with Z'^ = Z'^). 

Proof. See Appendix D. □ 

6.4 Lemma. Let v eW^, w eW. Then 

{v.dw)^0; (6.76) 
combined with Proposition 6.3, this gives v»dw E N^~^. 

Proof. For s = 1, ...,d, integration by parts and the assumption = divv = drV^ 
give 

(v. dw') - I dx v'drw' = ^ / dx (drv'')w' = . □ 

^ ^ (27r)'' Jt^ {2tt)'' Jt'' 
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6.5 Proposition. The bilinear map 

X ^ m^-' , (/, g) ^ /. dg (6.77) 

admits an estimate 

\\f.dg\\n-l^Kr,\\f\\n\\g\\n (6.78) 

for all f,g& IHI"^, with a suitable constant = K^a £ (0, +00). One can take for 
the latter any constant such that 



sup lCn{k) ^ Kr, , (6.79) 
fcezg 

^"^'^ (2.)^ J-', (1 + \hmi + \k- h\^r ^ ^ 

(as in (6.61), with Z'^ = Z^). 

Proof. See Appendix D. □ 
To conclude, we report another integral identity frequently used in the sequel. 
6.6 Lemma. For any v ^W^, one has 

{v\v.dv)L2^Q (6.81) 
Proof. We have (with q,s,re {1, d}), 

{v\v,dv)L2 = [ dxSgsv\v''drv') . (6.82) 

Prom here we infer, integrating by parts, 

{v\v»dv)L2 = — dx Sqsdr{v'^V^)v^ 

= -/" dxSgsidrV'iyv' - [ dxSgsV%dX)v' ^ - [ dxSgsidrv'^yv' (6.83) 

JTi JT^d J'Yd 

since div^; = 0; by comparison with (6.82) we obtain 

{v\v»dv)L2 — —{v\v»dv)L2 , (6.84) 
whence the thesis (6.81). □ 
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6.7 Remcirk. The inequality (6.73) (or (6.78)) is known from the hterature, in 
this form or in some variant: see, for example, one of the standard references on NS 

equations cited in the Introduction. To our knowledge, the novelty of Proposition 6.3 
(or 6.5) with respect to the already published material is the rule (6.74) (or (6.79)) 
to determine that can be used with Lemmas 6.1, 6.2 to provide a numerical 
value for this constant. Examples of this computation appear in Section 10 and 
Appendix H. 

The method employed in Appendix D to prove Propositions 6.3, 6.5 is very similar 
to one employed in [13] to estimate the product of two scalar functions in //"(R"^). 
In that paper, already mentioned in the Introduction, we have given a rule sim- 
ilar to (6.74) (or (6.79)) to find a constant Cnd = C'n such that ||pg||/fn(R<i) ^ 
C'n|bl|i?"(Rd)||?||i?"(Rd); furthermore, using convenient trial functions p,q we have 
shown this constant to be very close to the smallest one fulfilling the inequality. 
Due to the similarities with [13], the constant Kn provided by (6.74) (or (6.79)) is 
hopingly close to the smallest one for the inequality (6.73) (or (6.78)). 

Functions on T*^, depending on time. Suppose we have a function 

p : [to, T) C R ^ V or V , p{t) (6.85) 

where V ^ V stand for some spaces of R or R'' valued functions on T*^. At each 
"time" t G [to, T), this gives a function x e T'' i-^ p{t){x)] of course, we will use the 
more common notation 

p{x,t) ■.= p{t){x) . (6.86) 

The incompressible NS equations, in the Leray formulation. 

Let us recall that C°'^ indicates the locally Lipschitz maps. 

6.8 Definition. The incompressible NS Cauchy problem with initial datum vq G 
M^'^^ and forcing term rj G C'^'-'^([0, +00), H""^), in the Leray formulation, is the 
following. 

Find V G C([0, T),W^+^) n C^([0, T), H^-^), such that 
v{t) = Ai/(t) - Z{v{t). dv{t)) + r]{t) for t G [0, T) , = (6.87) 

(for some T G (0, +00] 

6.9 Remarks, (i) The requirement v G C([0, T), H""'"^) ensures by itself that the 
right hand side of the above differential equation is in C([0, T), EI^"^). 

(ii) The differential equation in (6.87) can be interpreted as the usual NS equa- 
tion for an incompressible fluid, in a convenient adimensional formulation. For each 
t G [0,T), v{t) : X G T*^ ^ v{x^t) is the velocity field of the fluid at time t; 
r]{t) : X >—>■ r][x, t) is the Leray projection of the density of external forces. Of course. 
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z/(t) is taken in the divergence free space EI^+^ to fulfill the condition of incompress- 
ibility; the pressure gradient does not appear in (6.87), having been eliminated by 
application of £. For completeness, all these facts are surveyed in Appendix E. 

The statement that follows refers to the time evolution of the functions t ^-^ {i^{t)) — 
(27r)~''/^d i'{x,t)dx and t ^ (l/2)||i/(t)|||2- Up to dimensional factors, the first one 
gives the mean value of the velocity field or, equivalently, the total momentum; the 
second one gives the total kinetic energy. 

6.10 Proposition. Suppose v fulfills (6.87) on an interval [0,T). Then, for all t 
in this interval, we have the following relations: 

(i) Balance of momentum: 

j^{v{t)) = m) . (6.88) 

(a) Balance of energy: 

^^ii^wiii^ = {m^^L^ + {mm)L^ ^ ^wi^w)^^ • (6.89) 

Proof, (i) Let us observe that 

= (z>) = (A^) - (£(^. dv)) + (r;) ; (6.90) 

the means of Ai/ and £(i/» dv) vanish due to Eqs. (6.12) (6.43) (6.76), so we get the 

thesis (6.88). 

(ii) Let us write 

\jtMl^ = (^l^)i'^ = - 9iy))L2 + {u\ri)L2 ; (6.91) 

on the other hand by the symmetry of the equality £z/ = u and Lemma 6.6, 

(z/|£(z/.(9z/))i2 = (£z/|z/.(9z/)i2 = {u\u,du)L2 = ; (6.92) 

these facts yield the equality in (6.89). The subsequent inequality in (6.89) follows 
via elementary integration by parts: 

{iy\Au)L2 = Srs I V/^v' = -5rs [ dv\dv' ^ . □ (6.93) 

Reducing the NS equations to the case of a zero mean velocity field. Let 

us recall the notation EI",j for the space of divergence free, zero mean vector fields 
(see Eq. (6.51) and subsequent comments); we regard this as a Hilbert space, with 
the inner product ( | )« and the norm || inherited from W^. 
The purpose of this paragraph is to show that the general Cauchy problem (6.87) 
can be reduced to a Cauchy problem for zero mean vector fields; let us define the 
latter precisely. 
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6.11 Definition. Let fo e H^+i and ^ e C^'^{[to,+oo),W^-^). The incompress- 
ible, zero mean NS Cauchy problem with initial datum vo and forcing term ^ is the 
following. 

Findipe C([0,T),e2+i) nCi([0,T),e2-i) such that 
ip{t) = A^{t) - £{^{t). difit)) + m for t e [0, T) , (^(0) = fo (6.94) 
(for some T e (0, +oo]). 

Let us connect this problem with the previous one (6.87), for given vq and rj. To this 
purpose, we need a bit more regularity on the forcing r^; to be precise, we assume 

voem^-^^ r]ec'''\[o,+oo),mi-^)nc{[o,+oo),w^) . (6.95) 

Let us define from vq and 77 the following objects: 

mo (vq) e R'^ ; me C([0, +00), R'^) , t ^ m{t) := mo + / ds {ri{s)) ; (6.96) 

Jo 

heC\[0,+oo),R'^) , t^h(t):^ [ dsm(s) ; (6.97) 

Jo 

fo:=Vo-moemi-^' ; (6.98) 

e : [0, T) ^ , t^ at) such that ^x, t) rj{x + h{t),t) - {rj{t)) . (6.99) 

The statement { (t) e H^^ for each t is evident from the definition. In Appendix F 
we will prove that 

eeC°'i([0,+oo),H^-^) . (6.100) 

6.12 Proposition. Let us consider: 

(i) the Cauchy prolem (6.87), with any datum Vq and forcing rj as in (6.95); 
(a) the above definitions of mo,m,h, fo,C^ and the Cauchy prolem (6.94). 
A function v of domain [0,T) fulfills (6.81) if and only if there is a function (f on 
[0,r) fulfilling (6.94) such that, for all x e and t E [0,T), 

iy{x, t) = m{t) + ip{x - h(t),t) . (6.101) 

Proof. Step 1. We suppose (6.94) with the above datum /o to have a solution (p on 
[0,T); we define v as in (6.101) and prove that it solves problem (6.87). It is clear 
that u is in the functional space prescribed by (6.87), and that the following holds: 

v{x, 0) = mo + <p{x, 0) = mo + foix) = vq{x) ; (6.102) 

z>(x, t) = {r]{t)) + <p{x - h{t),t) - h(t). d^{x - h{t),t) (6.103) 
= {r](t)) + Lp{x - h{t), t) - m{t), dcp{x - h{t)t, t) ; 
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Au{x, t) = Aip{x - h{t),t) ; (6.104) 

{u,du){x,t) = {m{t).d<^){x - h{t),t) + {^.d<^){x - h{t),t) . (6.105) 

The Leray projection £ commutes with space translations, so the last equation 
implies 

Z{iy,du){x,t) ^ £{m{t).d(p){x - h{t),t) + £{cp.dcp){x - h{t),t) . (6.106) 
To conclude, we note that 

£(m(i). dip) = m{t). dcp ; (6.107) 

to prove this, it suffices to check that m{t)»d(p is divergence free. In fact, 

dw{m{t).dip) = ds{m'{t)dr(p') = m''{t)dr{ds(p') = m"(t)a^(div(^) = , (6.108) 

since cp is divergence free at all times. Prom Eq. (6.102) we see that u fulfills the 
initial condition in (6.87); from Eqs. (6.103) (6.104) (6.106) (6.107) and (6.94), we 
see that u fulfills the evolution equation in (6.87). 

Step 2. Let us consider function u on [0,T), fulfilling (6.87); we will prove the 
existence of a function ip on [0, T) fulfilling (6.94), such that u and ip are related by 
(6.101). To this purpose, let us define a function cp by 

cp{x, t) := iy{x + h{t),t) - m(t) ; (6.109) 

then, at each time t e [0, T), 

{'p{t)) = {u{t)) - m{t) = (6.110) 

on account of Eq. (6.88) for (i/) and of the definition (6.96) of m. Besides having zero 

mean, (/? belongs to the function spaces prescribed by (6.87) due to the properties 
of u. Now, computations very similar to the ones of Step 1 prove that (p fulfills the 
initial condition and the evolution equation in (6.94). □ 

6.13 Remark. Let us regard Eqs. (6.98) (6.99) as defining a transformation 

r : W^+' X C°'i([0, +oo),Ul-') ^ X ^^'^([O, +oo),Ul;') , (6.111) 

(^^0,??) ^ {fo,0 = ^(vo,r]) . 
The map T is onto, due to the trivial equality (/o, ■C) = ^(/o, 0- 
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7 The NS equations in the general framework for 
evolution equations with quadratic nonlinear- 
ity. 

Basic notations. The zero mean version (6.94) is the final form for the NS Cauchy 
problem, to which we stick from now on. Let us recall that d e {2, 3,...}; throughout 
the section, we fix 

n>^, eeCO'i([0,+oo),H^-i) (7.1) 

(the function ^ is regarded to be given by itself, independently of any function rj 
as in the previous section). Our aim is to apply the formalism of Section 5 to the 
Cauchy problem (6.94) (and to the equivalent Volterra problem); in this case 

F± = (F±, II ||±) - (H^f , II |U±0 , F = (F, II II) := (H^„, || |U) , (7.2) 

A:=A: Mlt' ^ H^-i , f ^ Af ; (7.3) 
'J>:mi,xUl,-.W^-\ Tif,g):=-£{f.dg); ^ as in (7.1) . (7.4) 
Prom CP and ^, we construct the function 

V:W^,x [0, +oo) ^ e^-i , V(f, t) := 9{f, f) Hit) = -£(/. df) +m , (7.5) 

which appears in (6.94). 

For subsequent reference, we record the Fourier representations 

ii/iin±i = lY.i^+m-^w, ii/iin = 

(A/), = -\k\'fk . (7.6) 

Verification of properties (P1)-(P5'). The above set (F_,F,F+,^), completed 
with suitable semigroup estimators u, u_, has the properties (P1)-(P5) prescribed 
in Section 2, and (P4')(P5') of Section 5. We will indicate which parts of the proof 
are obvious, and give details on the nontrivial parts. 

7.1 Proposition. F=W^^, = W^f^ andA = A fulfill conditions (P1)-(P2). 
Proof. Everything follows easily from the Fourier representations. □ 

7.2 Proposition, (i) A generates a strongly continuous semigroup on H"^^, given 
by 

e*^/ = Yl e"*"' VfeSfc , forfe t e [0, +oo) . (7.7) 
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So, (P3) holds. 

(a) For f e and t e [0, +00) one has 

e*VeH^o, \\e'^f\\n^u{t)\\f%, u{t):=e-'; (7.8) 

the function {f,t) 1— >• e^^f gives a strongly continuous semigroup on H^q. 
(Hi) For f G IHI^i7^ '^'^^ ^ ^ (0; +00) one /ia5 

e*^/eH^, , ||e*^/|U^w_(t)||/|U_i , (7.9) 

{e^* 1 
^= /or < i ^ - , 
v/2it ^ 4 (7.10) 

V2 fort>-; 

note that u_{t), iJ.-{t) = 0{l/^/i) for t ^ 0'^. The function (/, t) ^ e*^/ is contin- 
uous from H^V^ X (0, +00) to H^^. 

With u_ as in (7.10), the function 1 1— >• U{t) := jQdsu_{s) is given by 

MIM) ^ 1 P-"' 



+ V2{e-^/^ - e"*) for - <t^ +00 



x/2 ' ^ 4 

7(i) := [ ds-^ forO^t^]. (7.12) 



/i •' 4 

(In particular, U{+oo) = 7(1/4)7^2 + ^6-^4 e (1.872, 1.873);. 
(^zf^ VFzt/i II _ as in (7.10), one has 

ft 



sup 

tG[0,+oo) 



[ ds ii_(t - s)e-' ^ V2 . (7.13) 
Jo 



(v) In conclusion, (P4)(P4') CLnd (P5)(P5') hold with 

B^l, a^\, N^V2. (7.14) 
Proof, (i) This follows basically from Eq. (7.6) for A. 

(ii)(iii) We only give details on the derivation of Eqs. (7.8)-(7.11). Let / G H^g, 
t e [0,+oo). Then, 

5^ (1 + \km{e'^f),\' = J] (1 + \k\'re-''\'\'\f,\' (7.15) 

^e-''Y.(^+\kmfk\\ 

fcGZg 

since |A;| ^ 1 for A; e Z^; this yields Eq. (7.8). Now, let / e H^"^ and t e (0,+(X)); 
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then, 

^ (1 + \k\y\{e''^ f),\' = E (1 + l^r)e-'*"l'(l + \k?T-'\fk? (7.16) 

^ ( sup u,{d))(Y,{^ + \k\'r-'\h?) , 



■i?e[i,+oo) 

Ut{^) := (1 +^?)e 







and an elementary computation gives 

{/ 1 \ e^* 1 

UJ 1) = — for < t ^ - , 
V2t J 2et 4 ' (7.17) 

Ut{l) = 2e-2* for i > - . 

Prom (7.17) we infer Eq. (7.9) with u_{t) := ^sup^g[i +(^) Ut{T&), i.e., 

e* 1 
for < i ^ 



\/2e-* for i > - ; 

this definition of u_{t) agrees with Eq. (7.10), and (7.11) follows trivially. 

(iv) See Appendix G. 

(v) Obvious consequence of items (i-iv). □ 

Analysis of IP and V. We turn the attention to the functions in Eqs. (7.4) (7.5). 

7.3 Proposition, (i) 7 is a bilinear map and admits the estimate (implying con- 
tinuity) 

\\nL9)\\n-i^Kjf\\Jg\U (7.19) 

for all f,g E W^q, with Kn = Knd any constant fulfilling (6.78) (so, condition 
(Ql) holds for this map). 

(a) As a consequence of (i), V fulfills the Lipschitz condition (P6). Furthermore, 
for each function (j) e C([0, T), H^^), the growth ofV from (f) admits the estimate 

\\V{f, t) - V{(t>{t)Mn-l ^ Lit, 11/ - m\n) (7.20) 

for t e [0, T) and f e M^,, where 
C:[0,+oo)x[0,r)^[0,+oo), (r, t) ^ C(r, t):=2X„||0(t) II r + X„r^ (7.21) 
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Proof, (i) The bilinearity is obvious, the estimate follows from Eqs. (6.49) for £ 
and (6.78) for the map (/, g) i— > /• dg. 

(ii) Use Corollaries 5.3, 5.4 on quadratic maps. □ 

The function H„_i. Prom here to the end of the paper, we denote in this way any 
function in C([0, +oo), [0, +oo)) such that 

5„_i e C{[0, +oo), [0, +00)) nondecreasing , (7.22) 

||e(t)||n-i^5n-i(^) forte[0,+oo) 
(e.g., En-i{t) := sup,g[o,4] ||C(s)||„_i). 
Cauchy and Volterra problems. 

7.4 Definition. For any fo e H^^^; CPn(/o) is the Cauchy problem (6.94), i-^-, 

Find ip G C([0,T),H;5+^) nC^([0,T),e^-i) such that 

if{t) = Aif{t) + V{^{t),t) for all t e [0, T) , (^(0) = /o . (7.23) 
For any fo e IHI"^, VP„(/o) is the Volterra problem 

Find if e C([0, T), M^J such that 

ip(t) = e*^/o + [ ds e^^-'^^V{ip{s), s) for all t e [0, T) . (7.24) 
Jo 

7.5 Remarks, (i) If fo e H^^^ VPn(/o) is equivalent to CPn(/o) by the reflexivity 
of the Hilbert space IHI^~^ (see once more Proposition 2.4). 

(ii) Por any fo G H";,, uniqueness and local existence are granted for VP„(/o) (Propo- 
sitions 2.5 and 3.10). 

8 Results for the NS equations arising from the 
previous framework. 

We keep the assumption (7.1) and all notations of the previous section; furthermore, 
we fix an initial datum 

fo e H^o ■ (8-1) 

The analysis of the previous section allows us to identify VPn(/o)with a Volterra 
problem of the general type discussed in Section 5, with semigroup estimators u, u_ 
of the form considered therein and a quadratic nonlinearity V. Due to Propositions 
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7.2, 7.3, the constant K and the functions u, u_, iJ,_, U, S_ of Section 5 can be taken 
as follows: 

X = a constant fulfilling (6.78) , u{t) := e"* , u_{t) = /i_(t)e-* , (8.2) 

as in (7.10) , W as in (7.11) , S_ = as in (7.22) . 
Hereafter, we rephrase Propositions 5.5 and 5.6 with the above specifications. 

8.1 Proposition. Let us consider for VP„(/o) an approximate solution (p^^ e 
C([0, T), H^q), where T e (0, +oo]. Suppose there are functions £n, T^^, TZn £ C([0, T), 

[0, +oo)) such that (i)-(iii) hold: 

(i) (pap has the integral error estimate 

\\E{p^.M\\n<Sr,{t) forte[0,T) ; (8.3) 

(a) one has 

yMWn^V^it) forte[0,T) ; (8.4) 
(Hi) TZn solves the control inequality 

£n{t) + Kr, [ dsu_(t - s)(2P„(s)7^„(s) +7^2(s)) ^ 7^„(^) forte [0,r) . (8.5) 
Jo 

Then, (a) and (b) hold: 

(a) VP(/o) has a solution if : [0,T) W^^; 

(b) one has 

y{t) - ^MWn < nn{t) forte [o,r) . (8.6) 

8.2 Proposition. Let us consider for VP„(/o) an approximate solution ipap ^ 
C([0, T), H^q), where T e (0,+oo]. Suppose there are functions Sn,T^n £ (^([O,?"], 
[0,+oo)) such that (i)-(iii) hold: 

(i) Sn is nondecreasing, and binds the integral error as in (8.3); 
(a) Vn is nondecreasing and binds (p^p as in (8.4); 
(Hi) one has 

2^ KnU{T)8n{T) + 2KJA{T)Vn{T) ^ 1 . (8.7) 
Then VP„(/o) has a solution ip : [0, T) EI",j and, for all t G [0, T), 

Mt)-Lp^M\n^nn{t) , (8.8) 

r l-2KM{t)Vn{t)-^{l-2KM{t)Vr,{t)y^KJA{t)8n{t) ^ ,^ ^, 

Un(0) ^/^ = 0; 

the above prescription gives a well defined, nondecreasing function TZn G C{[Q,T), 
[0,+oo)). 
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In Section 5, from Proposition 5.6 we have inferred Proposition 5.9, corresponding 
to the approximate solution := 0; in the present situation, this reads as follows. 

8.3 Proposition. Let 

3'n{t):^\\fo\\n + ^n-l{t)Uit) . (8.9) 

Suppose T G [0, H-cxo], and 

AKM{T)-Jn{T) ^ 1 . (8.10) 
Then VP„(/o) has a solution ip : [0, T) — > and, for all t G [0, T), 

Mt)\\n ^ Mt) X{AKJAmn{t)) (8.11) 

(where, as in (5.31): X{z) := py- — for z G (0, 1], A'(O) := I). 

The other results of Section 5 were about global existence and exponential decay, 
under specific assumption. In the present framework the constants S, o", N of the 
cited section are given by Eq. (7.14); this allows to rephrase Proposition 5.10 in this 
way. 

8.4 Proposition. Let us consider for VP„(/o) an approximate solution (p^^ G 
C([0, +oo), H"^). Suppose there are constants En,Dn G [0, +oo) such that: 

(i) admits the integral error estimate 

\\E{^a,){t)\\n ^ Enc-' for t G [0, +oo) ; (8.12) 
(a) for all t as above, 

yaMWn^D^e-' ; (8.13) 

(Hi) one has 



2 kpl^KnE^ + 2V2Kr,Dr, ^ 1 . (8.14) 
Then VP„(/o) has a global solution ip : [0, +oo) M^^ and, for all t G [0, +oo); 

Mt) - v^,{t)\\n ^ Rne-' , (8.15) 
_ 1 - 2V2KnDn - ^J{l- 2V2Kr,Dnf - AV2KnEn 

The apphcations of Proposition 5.10 considered in Section 5 were based on the as- 
sumption of exponential decay for the external forcing, that in the present framework 
must be formulated in this way: 
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(Q3)ji There is a constant J„_i e [0, +00) such that 

lleWlln-i ^ Jn-ie-^' for all t e [0, +00) . (8.16) 

The above mentioned applications in Section 5 were Proposition 5.12 (corresponding 
to the choice (p^^ := 0) and Proposition 5.14 (with ip^^ the ^-flow approximate 
solution). These can be restated, respectively, in the following way: 

8.5 Proposition. Assume ( Q3)n, and define 

i^n:=||/o||n + V2J„_i; (8.17) 

furthermore, assume 

4V2i^,F„ ^ 1 . (8.18) 

Then VP„(/o) has a global solution </? : [0, +00) — > H^^ and, for all t e [0, +00), 

Mt)\\n < FnX{AV2KnFr,)e-' (8.19) 

(with X as in (5.31)). 

8.6 Proposition. Define 

(/.„,eC([0,+oo),H^J , </P„,(t):=e*^/o+ f ds e^'-^^^^^s) . (8.20) 

Jo 

Furthermore, let us keep the assumptions and definitions (Q3)n (8.17) (8.18). Then 
the global solution tp : [0, +00) — >• H^g 0/ VP„(/o) is such that, for all t e [0, +00), 

Mt) - ^MWn < V2Kr,F^ X{AV2KnFn) e"* (8.21) 

(where, as in (5.60): X{z) := ^ ~ ^^^(^2^/^ ~ ^ -^^^ ^ ^ ^ 

8.7 Remarks, (i) Condition (8.10) can be fulfilled with cither /o, ^ small, T large 
or /o, ^ large, T small. Condition (8.18) is fulfilled if /o and ^ are sufficiently small. 

(ii) As a special case, suppose the external forcing ^ to be identically zero; then we 
can take S„_i = and J„_i = 0. Eqs. (8.10) and (8.18) ensure global existence 
if the datum fulfills the conditions 47^'„W(+oo) ||/o||„ ^ 1 and 4\/2A'„||/o||„ ^ 1, 
respectively. The less restrictive condition on /q is the second one, since \/2 < 
U{+oo). 

(iii) Let us return to Proposition 8.1; this states, amongst else, that the solution (/? 
of VP„(/o) exists on [0, T) if the inequality (8.5) has a solution TZ : [0, T) [0, +00). 
Let us compare this statement with a result presented in the recent work [2], that 
we rephrase here in our notations. 
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Let us consider the (incompressible, zero mean) NS equations with external forcing ^ 
and initial datum /o; when we refer to [2] a solution of this Cauchy problem means a 
strong solution, as defined therein. Now suppose </7„p to be an approximate solution 
on an interval [0, T) , and 

lkap(0)-/o||n-i ^ , \\ip^,it)-Aip^^{t)-Viip^M,s)\\n-i < en-i(i) , (8.22) 

||<^.p(t)IU-l ^ V,,_i{t) , ||v^„,(t)||„ ^ Vr,{t) 

for all t G [0,T), for suitable estimators ^ 0, e„_i, "D^ : [0,T) — > [0, +00). 

According to [2] (page 065204-10), the NS Cauchy problem has a solution ip : 
[0,T)^H2„ if 



rT ^ -Cn / ds (D„_i(sO + D„(sO) 

5n-i + / dsen-i{s) < 7777; e ^0 (8.23) 

Jo ^n-t 

(with C„ > a constant not computed explicitly, whose role is analogous to the one 
of Kn); (8.23) is an inequality involving only the approximate solution, and plays a 
role similar to our (8.5) to grant the existence of an exact solution on [0, T). 
Seemingly, Eq.(8.23) is not suited to obtain results of global existence for the exact 
solution. To explain this statement, suppose ip^p and its estimators to be defined 
on [0, -|-cx)), with 7^ or e„_i non identically zero; then (8.23) surely fails for 
large T, even in the most favourable situation where all integrals therein converge 
for T — > +00. In fact. 



T ^ +00 r+°° 
l.h.s. of (8.23) — > Sn-i+ dse„_i(s) e (0,+oo], (8.24) 

Jo 



1 +00 
< r.h.s. of (8.23) ^ — — — > . 

CnT 

On the contrary, our control inequality (8.5) can be used in certain cases to derive 
the existence of (/? (and bind its distance from Lp^p) up to T = +00; some applications 
of this type have appeared in the present section, further examples will be given in 
the next one on the Galerkin approximations (^^). 



^'^To conclude this remark we wish to point out that, under special assumptions, some global 
existence results could perhaps be derived from the approach of [2], with a different analysis of 
the differential inequalities proposed by the authors to infer Eq. (8.23). A discussion of this point, 
and of other interesting features of [2] , would occupy too much space here. 
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9 Galerkin approximate solutions of the NS equa- 
tions. 

Throughout this section, we consider a set G with the following features: 

G C , G finite , k e G ^ -k e G . (9.1) 

Hereafter we write -< >-keG ior the linear subspace of D spanned by the functions 
efc for k & G. 

Galerkin subspaces and projections. We define them as follows. 

9.1 Definition. The Galerkin subspace and projection corresponding to G are 
H^o ■= ®Eo^ -< ^keG^ VkCk \ Vk e ,v]^^ v^k ,k»Vk^ for all k} . 

keG 

(9.2) 

: D'^„ ^ H^o , '"=J2 ^V''v:=Y^ ^kCk ■ (9-3) 

fcezg keG 

It is clear that 

c n D',„ , A(efjce^„; (9.4) 

H^o C mZ , ?^(H^J = for all m e R . (9.5) 
The following result will be useful in the sequel. 

9.2 Lemma. Let n,p e R, n ^ p and v e Hso- Then, 

m-v>\\n<T^, \G\:= inf yrr^. (9.6) 

\G\P kez^\G 
Proof. We have (1 - ^*^)w = Zlfcezg\G ^^^fc' implying 

kez^\G kez^\G ^ ' ' 

whence the thesis. □ 
Galerkin approximate solutions. Let 

eeC([o,+cx)),D'^o) , (9.8) 

(of course, in the sequel V{f,t) := —£{f»df) + ^{t) whenever this makes sense). 
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9.3 Definition. The Galerkin approximate solution of NS corresponding to G, 
with external forcing ^ and datum fo, is the maximal solution c/?^ = (p^ of the 
following Cauchy problem, in the finite dimensional space H^p.- 

Find i{P e C\[0, Tg), M.g) such that (9.9) 

(^«(i) = A(^^(i) + ^«P((^«(i)) for all t, (^«(0) = ^«/o. 

Of course, "maximal" means that [0,Tg) is the largest interval of existence. In 
certain cases, one can prove that Tq = +00 and derive estimates of (p'^ (we return 
on this in the sequel). In this section we use the functions 

U as in Eq. (7.11) , A" as in Eq. (5.31) . (9.10) 

Let us consider any real number m, and assume ^ G C°'^([0, +00), EI™~^). We will 
use the notation S^_i to indicate any function in C([0, +00), [0, +00)) fulfilling Eq. 
(7.22) with n replaced by m; in the sequel that equation will be referred to as 
(7.22)^. When necessary we will make the assumption (QS)^ of exponential decay 
for the external forcing; this is like (Q3)„ with n — > m, thus involving a constant 
Jm-i e [0,+oo). 

The forthcoming proposition gives estimates on the interval of existence of the 
Galerkin solution and on its norm || \\m, which are in fact independent of G. 

9.4 Proposition. Let m > d/2, ^ e C^'\[0, +oo),W^-^) and fo G H™; then, 

hold. 

(i) Define (similarly to (8.9)) 

3'm(t) := \\fo\\m + ^m-l{t)U{t) ; (9.11) 

furthermore, letT e (0, +00], and assume the inequality 

AKJA{T)'J^{T) ^ 1 . (9.12) 
Then the Galerkin solution pP with this datum exists on [0, T) and fulfills 

Wv^'mn.^'D^it) forte[0,T), (9.13) 

-Dmit) := 3'm{t)X{4K^U{t)3'm{t)) . (9.14) 
(a) Alternatively, assume (Q3)m; define (similarly to (8.17)) 

i^m := ||/o|k + V2Jm-i , (9.15) 

and suppose 

4:V2KmFm ^ 1 . (9.16) 
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Then the Galerkin solution ({P with this datum is global, and fulfills 

ym^n < D^e-' forte [0,+oo) , (9.17) 

D^:^F^X{aV2K^F^) (9.18) 

(the above equations will be referred to in the sequel as (9.11)m, (9.12)m, etc.). 

Proof. We refer to the framework of Section 5 on systems with quadratic nonlin- 
earities. In the present case (F, || ||) := (H^g, || ||^), (F;p, || ||;p := (H^g, || ||m=Fi) 
have three copies of the same finite dimensional space, but equipped with differ- 
ent, though equivalent, norms); the operator ^ is A f H^^, and the bilinear map 
is : X eg, eg^, (/,^) v-^ <p^T(/,^); the function ^ of Section 5 is 

<pG^ G C°'i([0, +oo),egJ; the initial datum is <p^/o G eg,. 

For the operator A |~ Hg, we use the same estimates given for A in Proposition 7.2, 
with n — > m; this justifies using the scheme of Section 5 with U as in (7.11) and 

To estimate we use the inequalities on T in Proposition 5.2 with n ^ m, 

and the obvious relation H^"^ ■ ^ || ■ ||m-i; this gives ||^*^!P(/, 5()||to_i ^ 

-^m||/||m||5'||nn and so the constant K of Section 5 is, in this case, K^- 
For the initial datum ^^/o and for ^'^^ we use the estimates 

WfoWm < ll/olU , (9-19) 

||^«eWIU-i ^ \\m\\m-i < ^m-i{t) or J^_i ; (9.20) 

of course, the bound via Em-i refers to case (i) and the bound via Jm is for case 
(ii). Applying to this framework Propositions 5.9 and 5.12 we get the statements in 
(i) and (ii), respectively. □ 

9.5 Remark. Global existence of ip'-^ could be proved under much weaker con- 
ditions than the ones in item (ii) of the above proposition. In fact, using for 
an energy balance relation similar to (6.89), one can derive global existence and 
boundedness of ||<^'^(i)||2;,2 when /q is arbitrary and the external forcing makes finite 
both integrals jQ°°dt \\C{t)\\L'2, j^^^dt \\i{t)\\\2 ■ see, e.g., [15]. However, the ener- 
getic approach does not allow to derive estimates of the specific type appearing in 
Proposition 9.4. 

The distance between the exact NS solution and the Galerkin approxi- 
mations. From here to the end of the paragraph, we fix two real numbers 

P>n>^; (9.21) 
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we also fix 

eeC°'^([0,+oo),H5-^), foemi, (9.22) 

and denote with (p*^ the Galerkin approximate solution with such forcing and datum, 
for any G as before. This will be compared with the solution (p of the NS equations 
with the same forcing and datum. 

9.6 Lemma. Let us regard as an approximate solution o/VP„(/o) ; then the 
following holds. 
(0) The integral error of ip^ is 



E{<p^){t)^-{1-^' 



G\ 



Jo 



(9.23) 



(i) Let us introduce the definitions or assumptions (9.11)p (9.12)p, for some T e 
(0, +oo] (implying existence of ^pP on [0,T)j. Then, for all t e [0,T) we have 



ML 

\Q\v-n ' 



yp{t) :=5-p(t) \^ KpU{t)'5,{t)X\^K^{t)'Jp{t)) 



(9.24) 



(9.25; 



The function is nondecreasing. 

(a) Alternatively, introduce the definitions or assumptions (Q3)p (9.15)p (9.16)p 
(implying that ^pP is global). Then, for all t e [0, +oo) we have 



\G\p~"' 

1 + V2KpFpX^{4:V2KpFp)^ 
Proof. Derivation of (9.23). By definition 



y •= F 
J. p . ± p 



(9.26) 
(9.27) 



E{^^){t) = ^«(t) - e*^/o - rdse(*-)^P(^^(5)) ; (9.28) 

^0 

on the other hand, the Cauchy problem (9.9) defining ipP has the integral reformu- 
lation ^ 

(^^(i) = e*^^^/o + / dse(*-^)^q3^7'((^^(s)) ; (9.29) 
Jo 

inserting this into (9.28) we get 



E{^'^){t) = -e*^(l - ^^)/o - f dse^'-'^^{l - ^^)P(99^(s)) . 

Jo 



(9.30) 
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Finally, the operator 1 — commutes with A and its semigroup (as made evident 
by the Fourier representations); so, 1 — can be factored out and we obtain the 

thesis (9.23). 

Some preliminaries to the proof of (i) and (ii). From (9.23), the estimates (9.6) on 
1 — and (7.8)-(7.10) on the semigroup of A we get 



\Q\p-n 



Jo 
Jo 



(9.31) 



Ip-i 



Proof of (i). From P(<^*^(s)) = '?{ip'^{s),ip'-'{s)) + ^(s) we infer the following, for 
s e (0,t): 

\\Vi<f%))\\,., ^ K,\\<f^is)\\l+ms)\\p-i (9.32) 

^ K.Vlis) + E,_,{s) ^ K,Vl{t) + Sp_i(t) 

(in the above, we have used (9.13)p (7.22)p and the relation T>p(ys) ^ T>p{t)). 
We insert the result (9.32) into Eq. (9.31); in this way we are left with an integral 
Jg*dse~(*~*)//_(i — s) — ds fj,_{s) ^ l{{t). From this bound and e~* ^ 1 we 
obtain 



\Q\p-n 



\G\P- 



||/o||p + W(i)(XpDj(t) + ViW) ] (9-33) 



where the last passage follows from definition (9.11)^; now, explicitating Vpit) we 
get the thesis (9.24) (9.25). Finally, yp is nondecreasing because J^p, U and X are 
so. 

Proof of (ii). In this case, from the inequality ||P((^'^(s))||p_i ^ i^p||(^'^(s)||p +||^(s)||p_i 
we infer, by means of Eqs. (9.17)p and (Q3)p , 



\\V{<p''{s))\\p.,^KpDle-^^ + Jp^,e 



(9.34) 



Inserting this result into (9.31) we are left with a term e ^ J^dse " li-{t — s), which 
is bounded by \/2e~* due to (7.13); the conclusion is 



\Q\p-n 



\G\v 



\\h\\p + V2{KpDl + Jp_,) 
Fp + V2KpDl 



(9.35) 
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(the last equality following from (9.15)p). Now, explicitating Dp we get the thesis 
(9.26) (9.27). □ 

The following proposition contains the main result of the section. 

9.7 Proposition, (i) Let T G (0, +oo]; make the assumptions and definitions 
(9.11)n (9.12)n and (9.11)p (9.12)p (implying the existence of^p''' on [0,T)). Finally, 
with T>n andyp defined by (9.14)n (9.25), assume 



+ '^KM{T)V^{T) ^ 1 . (9.36) 
Then VP„(/o) has a solution (/? of domain [0, T) and, for all t in this interval, 

M^)-^^WIIn^^^^^#, (9-37) 

The function t ^ Wnp |G|(^) ^■s nondecreasing; a rough, \G\-independent bound for it 
is 

'^"°lW S-2wjp„(T) <8-^8) 

for allt e [0,T). 

(a) Alternatively, make the assumptions and definitions (Q3)n (9.15)n (9.16)n and 
(Q3)p (9.15)p (9.16)p (implying global existence of ip^). Finally, with Dn and Yp 
defined by (9.18)^ (9.21), assume 



2 + 2^/2i^„i?n ^ 1 . (9.40) 

Then VPnifo) has a solution ip of domain [0, +oo) and, for all t in this interval, 

Mt)-(^«(t)|U^|^e-\ (9.41) 

Wnp :^ ^ X ( ] . (9.42) 

l-2V2KnDn \(l-2V2KnDnr\G\P--J ^ ' 

The above constant has the rough, \G\-independent bound 
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Proof, (i) A simple application of Proposition 8.2, with 

= P„ as in (9.14) ; Sn{t) = . (9.44) 

The condition (8.7) in the cited proposition takes the form (9.36). The proposition 
ensures that (p is defined on [0,T), and gives the estimate 

Mt)-cp''{t)\\n^nn{t) , (9.45) 

involving the nondecreasing, continuous function 

^ ... ._ 1 - 2KM{t)Vn{t) - V(l - 2KM{t)Vn{t)Y - 4.KM{t)yp{t)/\G v-- 



2X„W(i) 

^"(0) ^= 1^ ■ (9-46) 



We note that we can write 



(1 - 2ir„W(t)P„(t))|G|f-« V(l - 2ii'„W(t)P„(t))2|G'|f 
this yields the thesis (9.37) (9.38). 

The fact that VVrep|G| is a nondecreasing function of time is apparent from its defini- 
tion. The bound (9.39) for it follows from the nondecreasing nature of the function 
t ^ yp{t)/{l - 2KJJ{t)Vn{t)) and from the inequality X{z) ^ 2 for all z e [0, 1]). 
(ii) A simple application of Proposition 8.4, with 

(^., = <^^, D„ as in (9.18), K = ■ (9.48) 

The condition (8.14) in the cited proposition takes the form (9.40). The proposition 
ensures that (f is defined on [0, +oo), and gives the estimate 

\\ipit)-ip''it)\\n^Rne-' , (9.49) 



1 - 2V2K^D^ - J{1- 2V2i^„L»„)2 - AV2K^Yp/\G\P-- 
■ (9-50) 

We note that we can write 

Rn = X ^ ^ , 9.51 

{1-2V2KM\G\P'- \{1-2V2KM''\G\P--J 

yielding the thesis (9.41) (9.42); the rough bound (9.43) follows again from the 
inequality X{z) ^2. □ 
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9.8 Remcirk. Of course, if p > n the previous proposition implies convergence of 
the Galerkin solution to the exact solution 99 of VP„(/o) . More precisely, with 
the assumptions in (i) we infer from (9.36) (9.39) that 



sup \Mt) - cp^{t)\\n = 0{—^) ^ for IGI ^ +00 ; (9.52) 



te[o,T) 



\Q\p-n 



sup e'Mt)-<p''{t)\U^O{—^)^0 forlGH+oo. (9.53) 



in case (ii), we infer from (9.41) (9.43) that 
sup e'Mt)-v''mn^O{-^^ 

10 Numerical examples. 

Given the necessary constants Kn, the datum norms and some bounds on the ex- 
ternal forcing, the framework of Sections 8 and 9 yields informations on the time of 
existence of the solution ip of VP„(/o) , and on its distance from an approximate 
solution. In the sequel we exemplify such estimates referring to Section 9, i.e., to 
the Galerkin approximations. 
Throughout the section, we take 

d = 3; n = 2, p = 4. (10.1) 

The constants K2 and K4 involved in calculations can be obtained from Lemmas 6.1 
6.2 and Proposition 6.3; a MATLAB computation illustrated in Appendix H yields 
the values 

K2 = 0.20 , K4 = 0.067 . (10.2) 

The other calculations mentioned hereafter have been performed using MATHE- 
MATICA. 

An application of Proposition 9.7, item (i). We suppose the external forcing 
has bounds (7.22)2 (7.22)4 with Si(t) = const. = Si and '^sit) = const. = S3 for 
all t e [0, +00). Conditions (9.12)2 and (9.12)4 are satisfied with T = -|-cxd if 

||/o||2 + 1.88Si< 0.667, II/0II4 + 1.88 S3 < 1.99 ; (10.3) 

under the above inequalities for /o and the forcing, the Galerkin solution (p'^ exists 
on [0, +00) for each G. As an example, conditions (10.3) are satisfied in the case 

II/0II2-O.I5 , ||/o||4 = 1.50. Si = 0.025, S3 = 0.25, (10.4) 

to which we stick hereafter. In the above case, condition (9.36) with n = 2, p = 4 
and T = +00 becomes 0.161 + 2.31/1^1 ^ 1, which is fulfilled if 

IGI ^ 2.76 . (10.5) 
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Assuming (10.4) (10.5) the solution of VP2(/o) is also global, and 

Mt) - (^«(i)||2 ^ ^^^^0^ < ^ for all G as above, t e [0, +oo) . (10.6) 

The numerical value of W24|g|(^) can be computed at will from definition (9.38); 
here we have used the rough bound VV24|g|(^) ^ 8.71, coming from (9.39). 

Another application of Proposition 9.7, item (i). We maintain the assump- 
tions Si(t) = const. = Si and S3(t) = const. = S3 for all t G [0, +cxd). We 
take 

||/o||2 = 0.20 , ||/o||4 = 2.00. Si, S3 as in (10.4). (10.7) 

Now conditions (10.3) are not fulfilled, indicating that (9.12)2 and (9.12)4 are not 
satisfied with T — +00. On the contrary, (9.12)2 and (9.12)4 are found to hold with 

T=1.51, (10.8) 

i.e., the Galerkin solution (p'~^ exists for any G on the time interval [0,1.51). To 
go on, we note that condition (9.36) with n — 2, p ~ A and T as above becomes 
0.163 + 2.41/|G'| ^ 1, which is fulfilled if 

IGI ^ 2.88 . (10.9) 

Under the assumption (10.9) the solution (f of VP2(/o) exists on the same interval, 
and 

Mt) - (^«(i)||2 ^ ^ ^ for all G as above, t e [0, 1.51) . (10.10) 

I I I I 

Again, we can compute the numerical value of W24|g|(^) from the definition (9.38); 
here we have used the bound W24|g|(^) ^ 11.1, coming from (9.39). 

An application of Proposition 9.7, item (ii). Let us recall that this case 
refers to exponentially decaying forcing. From the datum norms ||/o||m and the 
constants Jm-i in the forcing bounds, as in (9.15) we define the coefficients :— 
||/o||m + V^Jm-i for m = 2,4. Conditions (9.16) for m = 2,4 become, respectively, 
I.I4F2 ^ 1 and O.38OF4 ^ 1; these are fulfilled if 

F2 < 0.877 , F4 ^ 2.63 , (10.11) 

and in this case the Galerkin solution (p'-^ is global for each G. As an example, let 
us suppose 

F2 = 0.20 , F4 = 2.00 ; (10.12) 
then, condition (9.40) becomes 0.121 + 1.75/|G| ^ 1, which is fulfilled if 

IGI ^ 2.00 . (10.13) 



61 



With these assumptions the exact solution (/? of VP2(/o) is global, and 

Mt) - (^«(i)||2 ^ ' ^ for all G as above, t e [0, +00) . (10.14) 

\G\ \G\ 

The expression of 14^24 101 is provided by (9.42); here we have used the rough bound 
W^24|G| ^ 6.10, coming from (9.43). 



A Appendix. Proof of Lemma 2.6. 

First of all, we put 

Z := sup z{t) ; (A.l) 

t&[to,T] 

we continue in two steps. 
Step 1. For aJi A; e N, one has 

m g z (oneM. (A.2) 

To prove this, we write (A.2)fc for the above equation at order k, and proceed by 
recursion. Eq. (A.2)o is just the inequality z{t) ^ Z. Now, we suppose that {A.2)k 
holds and infer from it Eq. (A. 2)^+1. To this purpose, we substitute (A.2)fe into the 
basic inequality (2.22), which gives 

r(ka + l) i,, ((-5)1-" 

expressing the integral via the known identity (4.20), we get the thesis (A. 2)^+1. 
Step 2. z{t) — for all t e [to^r]. In Eq. (A.2), let us fix t and send k to oo; the 
right hand side of this inequality vanishes in this limit, yielding the thesis. □ 



B Appendix. A scheme to solve numerically the 
control inequality (5.14). 

Notations. In this Appendix we often write {0, ...,M} where M is an integer or 
+00. If M is a nonnegative integer this will mean, as usually, the set of integers 
0, 1, 2...., M. If M is a negative integer, we will intend {0, M} := 0. If M = +00, 
{0, M} will mean the set N of all natural numbers. 
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We often consider finite or infinite sequences of real numbers of the form {tm)me{o,...,M}'i 
if M = +00, we intend tj^ '■— ^^^m-^+oo tm whenever the hmit exists. 
The numerical scheme. Let us be given an approximate solution cp^^ e ^"([0, T'),F) 
of VP(/o) , where T' e (0, +00]; in the sequel we always intend 

£{t) := , V{t) := ||^.,(t)|| for t G [0,T') . (B.l) 

Hereafter we outline a numerically implementable algorithm to construct a solution 
TZ of the integral inequality (5.14) on some interval [0, T) C [0, T'); this solution TZ 
will be piecewise linear. 

In order to construct the algorithm, we choose a sequence of instants {tm)m=o,...,M', 
where M' is a positive integer or +00. We assume 

= to <ti <t2 < ... <tM' = T' . (B.2) 

Furthermore, we denote with £m, I'm, Hmk, ^mk, N^k some constants such that 

sup £{t) ^ Era , sup V{t) ^ ; (B.3) 

i G \pm 5^771+ 1 ) is \tm )im+ 1 ) 

' ^ rr ..... T'^^.. .^ ^-^fc / , 

mk 



sup I ds u_{t-s) ( — ^ i^mfc, sup f dsu_{t-s)-^ — ^^7, 

ftk+l 

sup / ds u_{t - s) ^ Nrnk for m e {1, ...,M' - 1}, A; e {0, ...,m - 1} ; 

sup dsu_{t-s) (-^ — ~r] ^ dsu_{t-s)--^ — 3^ ^ /m 

te[tm,tm+l) Jtm \tm+l J te[tm ,tm+l) J tm 

sup jds u_{t -s) ^ Nmm for m e {0, M' - 1} . (B.4) 
Finally, for m, k as above and all a, x e R, we define 

^mk{a, X) := {Hmk + Nmk - "^IrnkW + 2(/mfe - Hmk)aX + Hrnkx"^ (B.5) 

+2(A^mfe - Imk)T^ka + 2ImkT^kX . 

B.l Proposition. Suppose there is a finite or infinite sequence of nonnegative 
reals (7^m)me{o,...,M} (with 1 ^ M ^ M') such that 

m 

£:^ + xJ]$^fc(7^fc,7^fc+l) ^min(7^^,7^^+l) forme {0,....,M -1} . (B.6) 

fe=0 
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Let TZ e C([0, ^m), [0,+oo)) be the unique piecewise linear map with values TZ^ at 
the times tm, i-e., 

n{t) = TZm+C^m+l-nm) ^"^"1 fortE [tm,tm+l), m G {0, M - 1} . (B.7) 

Then, TZ solves the integral inequality (5.14) on [0,tM)- 

Proof. Let TZ be defined as above, and t in some subinterval [tm,tm+i) £ 
{0,...,M- 1}). Then 

l.h.s. of (5.14) ^Sm + K [ dsu_{t - s){2V{s) + 7^(s))7^(s) (B.8) 

JQ 

^£m + K\Y^ I \dsu_{t- s){2Vk + n{s))n{s) 



^S^ + K f"' + f ]dsu_{t-s) (2Vk + 7^fc + {TZk+i - 7^fc) / \ ] x 

Vfe^-^*fc -^tm/ V tk+l-tkj 

(— f \ 
tk+i-tkj 

the last passage following from the inequalities (B.4) and the definition (B.5) of ^mk- 
From here and from (B.6) we infer, for t in the same interval, 

l.h.s. of (5.14) ^ min(7^^,7^^+l) ^ TZ{t) = r.h.s. of (5.14) (B.9) 

In conclusion, (B.6) ensures TZ to fulfill (5.14) on [0,tM)- D 



B.2 Remcirks. (i) A sequence of constants {£k) fulfiUing the first inequality (B.3) 

is easily obtained if ip^^ G C([0, T'), F+) D C^([0, T'), F_), and there are suitable 
estimators for (the semigroup and) for the datum and differential errors d{ipap), 
e{ip^p). More precisely suppose that 

\\d{<pJ\\^S, (B.IO) 

and that, for m G {0, M' - 1}, 

sup U{t) ^ Um , sup ||e(</7ap(*))||- ^ Cm , (B-ll) 

(iim)m=o,...,M'-i and (6^)^=0,. ..,M'-i being sequences of nonnegative reals. 
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Prom Lemma 3.4 on the integral error we obtain, for t e [tm, tm+i), 

^ +/ \dsu_{t-s)\\e{^^M)\\-; (B.12) 

now, from (B.ll) and (B.4) we infer, for m G {0, 1, M' — 1}, 

m 

£{t)^£rn a t e [tm,tm+l) , £m ■= UmS + ^ Nmk^k ■ (B.13) 

k=0 

(In fact, one could extend this result to the case where ip^p is continuous from [to, T') 
to F+ and piecewise from [to, T') to F_: this typically occurs for the approximate 
solutions defined by finite difference schemes in time). 

(ii) For any m, Eq. (B.6) holds if and only if 

m 

either 7^^+l e [0, Um) , £m + Kj2 ^fc+i) ^ ^m+i , (B.14) 

k=0 

m 

or n^+i e [7^^, +oo) , £m + K^ ^mk{T^k: T^k+i) ^ T^m ; (B.15) 

fc=0 

note that each ^^nk in these inequalities is a quadratic polynomial. For m = 0, Eq.s 
(B.14) (B.15) define a problem for two unknowns {1Zq,1Zi); for m > 0, we can see 
(B.14) (B.15) as a problem to determine recursively TZm+i from TZq, ...,TZm- 
For each m, if the problem has solutions it seems convenient to choose for TZm+i 
the smallest admissible value. This criterion could be applied for m = as well, 
choosing among all solutions {TZo,TZi) the one with the smallest TZi. 

(iii) In practical computations, the determination of 7^„,+i from TZq, ...,7^m goes on 
until problem (B.14) or (B.15) have solutions. The iteration ends if, for some finite 
M, both the above inequalities for TZm+i have no solutions. In this case, we have a 
function TZ solving (5.14) on the interval [0,tM)- Alternatively, the iteration might 
go on indefinitely. 

(iv) The recursive scheme (B.6) has a typical feature of the iterative methods to 
solve integral equations or inequalities of the Volterra type: to find 7^m+i one must 
compute a "memory term" involving 7^0i---^m- The memory term depends non- 
trivially on m (through the coefficients Imk, etc.), so it must be fully redetermined 
at each step; this makes the computation more and more expensive while m grows. 
An exception to this framework occurs if the semigroup estimator t i— > u_(t) is (a 
constant x) an exponential, at least for t greater than some fixed time -(9; this is 
just the case of the NS equations, see the forthcoming Remark B.4 (ii). In this 
special situation, Eq. (B.6) can be rephrased as a pair of recursion relations for 
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two real sequences {TZm), {Sm)] at each step, computation of TZm+i and »S^+i does 
not involve the whole previous history, but only the values of Sm and of TZk for 
tm — < tk+1 ^ tm+1- The forthcoming Proposition explains all this. 

B.3 Proposition. Let us suppose there are i? ^ and A,B > such that 

u_{t) = Ae-^^ for t e {-d, +00) ; (B.16) 

furthermore, let us intend that k always means an integer in {0, M' — 1}. Then, 
i) a) hold, 
i) Let 

mG {0,...,M'} , tfe+i^t^-^?; (B.17) 



then, conditions (B.4) are fulfilled with 



B^tk+i - tk) 



2 



-'mfe ■ 1 ■ 7-)9/i I \ ' 

B\tk+i - tk) 

TV^fe := Ae-^'-Nk , := . (B.18) 

Consequently, for all real a, x one has 

^mk{a, x) = ^e-^*'"$fe(a, x) , (B.19) 

$A:(a, x) := {Hk + Nk- 2h)a^ + 2(4 - Hk)ax + i/fcX^ + 2(iVfc - 4)Pfca + 24Pfca; . 

a) Consider a sequence (7^m)me{o,...,M} of nonnegative reals. Then, (TZm) fulfills Eq. 
(B.6) if and only if there is a sequence of reals (5^)me{o,...,M-i} such that 

Sm+ M'^k,nk+i)^Sm+i forme{0,...,M-2} ; (B.20) 

{k\tm.-^<tk+l^tm+l-^} 

S^ + K J2 $mfc(7^fc,7^fe+l) + XAe-^*-5^ ^min(7^^,7^^+l) (B.21) 

{k\tm-'6<tk+l^tm.+ l} 

forme {0, ...,M- 1} . 

Hi) In particular, suppose the instants tm and 'd to he integer multiples of a basic 
spacing r > 0; 

t^ = mr for m e {0, M'} ; -d = Lt for some LeN . (B.22) 
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Then, for m e {0, M - 1}, 

{k\tm-^<h+i^tm+i}-^ {0,...,m} ^fm<L, ^^'^^^ 

{k\tm - ^ < t.+i ^ Wi - n = { " IZ<l[ (^-24) 

Proof, i) Let t G [tm,tm+i)- For s e [tk,tk+i) one has i — s > — ^fe+i ^ 
implying — s) = Ae~^^^~^^; so, 

f y = Ae-^* Tds e^^ f ^ ~ V = Ae'^Hk (B.25) 

Jtfc \tk+i-tkj Jtk \tk+i-tkj 

^ Ae"-^*'"i/fe = i/^fc as in (B.18) . 

In conclusion, defining Hmk as in (B.18) we fulfill the first inequality in (B.4) (note 
that k — 1 due to tk+i ^ tm)- Similarly, the other inequalities (B.4) are fulfilled 
with Imk, Nmk as in (B.18). 

Finally, inserting Eq. (B.18) into Eq. (B.5) for wc obtain the thesis (B.19). 

ii) Let us rephrase Eq. (B.6) in the case under examination. To this purpose, we 
reexpress the sum therein writing 

m 

k=0 {k\tm-'fi<tk+l^tm+l} {fe|tfc+l^tm— 1?} 

and then use Eq. (B.19) for the summands with tk+i ^ tm — i^; in this way, Eq. 
(B.6) becomes 

{k\tm-'&<tk+l^tm+l} {kUk+l^tm-l^} 

^ min(7^™, Um+i) for m G {0, M - 1} . (B.27) 

Let us consider any sequence {'R-m)me{o,...,M} nonnegative reals. If {TZm) fulfills 
(B.27), define 

Sm-= 5Z '^^•(^'^'^^+1) formG {0,...,M-1} ; (B.28) 

{k\tk+l^tm-&} 

then (B.20) (B.21) follow immediately (in fact, with ^ replaced by = in (B.20)). 
Conversely, suppose there is a sequence of reals (iSm)mG{o,...,Af-i} fulfilhng Eq.s (B.20) 
and (B.21) with (TZ^); then "^fe+i) ^ ^5^, and it is easy to 

infer Eq. (B.27) for (7^^). 

iii) Obvious. □ 
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B.4 Remcirk. (i) Eq. (B.20) does not prescribe ^o- It is convenient to choose 
So '■— 0; with this position Eq. (B.21) with m = is a problem for two un- 
knowns (7^0) T^i), for which we could repeat the comments of Remark B.2 (ii). After 
these initial steps, we can use Eq.s (B.20) (B.21) as recursion relations to obtain 
iSi, 7^.2, ^2, 7^3, and so on. 

(ii) As anticipated. Proposition B.3 can be applied to the NS equations, in the 
framework of Section 7. Eq. (7.10) of the cited section gives the semigroup estimator 

u_{t) := e*/\/2ei for t ^ 1/4, and u_{t) := A/2e"* for all t > 1/4. So, the conditions 
of the previous proposition are fulfined with = 1/4, A = B = 1. These values 
must be sustituted into Eq. (B.18) for Hmk, Imk and Nmk giving, for example, 

AT^fe := \/2e-'"^Nk , Nk := e*'=+i - 6**= for tk+i ^tm- 1/4 . (B.29) 

Explicit expressions could be derived as well for H^^, and when tm — 
1/4 < tk+i ^ tm, using elementary bounds on u_ derived from the expression (7.10). 
However, this requires a tedious analysis of a number of cases, since the parameter 
t — s in Eq. (B.4) can be smaller or greater that 1/4; details will be given elsewhere, 
when we will treat systematically the approach outlined in this Appendix. 



C Appendix. Proof of Lemmas 6.1 and 6.2. 

We begin with two auxiliary Lemmas. 

C.l Lemma. Let us consider two radii p, pi such that 2\fd ^ p < pi ^ +00, and 
a nonincreasing function % e C([p — 2^/d, pi), [0, +00) ). Then, 

o d/2 fpi 

E ^(1^1) ^ WJm / + ^y~'x{t) , (C.l) 

Proof. Let us introduce the cubes 

/i + [o,i]'^ {heZ"^) (C.2) 

and the annulus 

A{p-Vd,pi + Vd):={qeR'^ \ p-Vd^\q\<Pi + Vd} ; (C.3) 
we claim that 

UfteZ'*,p^H<pi(/i+[0,l]')c A(p-v^,pi + ^) . (C.4) 

In fact, q e h + [0, 1^ implies \h\ — Vd ^ |g| ^ \h\ + y/d; now, if p ^ \h\ < pi we 
conclude p — \/d ^ \q\ < pi + \/d. 
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The inclusion (C.4) implies 

E / dqx(\q\-Vd)^ f dqx(\q\-Vd). (C.5) 

hez'i,p^\h\<p, '^'^+[0'il' JAip-Vd,p,+Vd) 

On the other hand, for h as in the above sum and q e h + [0, l]"^, we have \q\ — Vd ^ 
\h\, whence x{\q\ ~ v^) ^ xd^Di "this implies 

/ dqx{\q\ -Vd)^ x{\h\) f dq = x{\h\) . (C.6) 

Jh+[0,l]<^ Jh+[0,1]^ 

Prom here and from (C.5) we obtain 

E [ dqx{\q\-Vd) . (C.7) 

hezd,p^\h\<pi J A{p-Vd,p^+Vd) 

The right hand side of (C.7) can be expressed in terms of the one-dimensional 
variable r — \q\; as well known, dq — (27r'^/^/r(d/2))r'^~^dr, so 

hez'i,p^\h\<pi ^ ' ' -^p-"^ 

now, a change of variables t — r — Vd gives the thesis (C.l). □ 
To go on, let us recall the convention (6.56) 

Z'^ := Z'^ or . 

C.2 Lemma, (i) Generalizing (6.59), let 

^ TTtW for^A^i^.^^). (C.9) 
^ ' h&Z<i,\h\<\ ^ ^ ^ ^ > 



Then, 



§.(A) ^S.+ 2^-%^l4(d/2) ^-^^^ /or.>0,A>2v^, (C.IO) 

s V ^ 

(271)'^ ^ (l+|/i|2)-' 

/ieZrf,|/i|<2\/d 



log(^) z/. = d/2. 
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For fixed u and A — > +00, this implies 



0(1) tfu>d/2, 
§,(A)=<( O(logA) ifu = d/2, (C.U) 
0(A'='-2^) tf0<u<d/2. 



(a) Let 



^^^^^^ '■= (i^^ ^ (TtW /«-^>^A^>0- (C-12) 

r/ien 

A§^(A) ^ 5§^(A) /or 1/ > d/2, A > 2\/d , (C.13) 
where (generalizing (6.60)) 

5§.(A) := ^^^-^ p . (C.14) 

Proof, (i) Let > 0, A > 2\fd. Dividing in two parts the sum defining §t/(A), we 
get 

hGZ<*,2%/d^|h|<A 

Now, we bind the sum using (C.l) with x(t) := 1/(1 + t^Y , p := 2\fd and p\ := A, 
yielding 

1 ^ [t + ^rf-^ 



On the other hand, one estabhshes by elementary means the inequality 

1 ^ 1 + d ^ ^ 



l + f" {t + Vdy 

(holding as an equality when t — 1/y/d). Inserting this into (C.16), we get 

(I + 1 

Sz.(A) ^ 5^ + -r-r-n^T;TTTT / dt p . (C.18) 

- 2<^-%W(d/2) Jo (t + ^)2-'^+i ^ ^ 

The last integral equals Fi,{X), so (C.IO) is proved. Having this result, the statement 
(C.ll) on the limit A — > +00 is obvious. 

(ii) Let u > d/2, A > 2Vd. To bind AS„(A) we use (C.l) with xit) := 1/(1 + f^Y, 
p :— X and pi :— +00, and subsequently employ the inequality (C.17). This gives 
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1 r+°° (f _L ./^)d-i 



^ ; / dt ■ 



and computing the last integral we justify Eqs. (C.13) (C.14). □ 

From here to the end of the Appendix we fix a real number n, fulfilling the relation 
(6.55) 

d 

here are the proofs of Lemmas 6.1, 6.2. 

Proof of Lemma 6.1. We take any A ^ 2-\/d; with the notations of the previous 
Lemma, we have 

S„ = S„(A)+A§„(A) . (C.20) 

Both terms in the right hand side are finite; the term A§„(A) has the upper bounds 
(C.13) and the obvious lower bound AS„(A) > 0. From these facts we infer the 
finiteness of S„, and the bounds (6.58) for it. □ 

Proof of Lemma 6.2. We consider any cutoff function A. : 7/ ^ [2-\/d, +oo). For 
k e 7/, we introduce the decomposition 

K,n{k)^%n{k) + A%n{k) ; (C.21) 

here %n{k) is defined by (6.66), and 

For the term 3C„(/i;), we furtherly introduce a decomposition 

%r.{k)^%'M + K{k) , (C.23) 

(27r)'^ ^ (l+|/,|2)n(l+|^_/,|2)n \^-^^) 

(the sum defining 3C"(A;) is meant to be zero if A(/c) ^ \k\l'l). In the sequel we 
analyse separately A3C„, and 3C^; we will frequently use the inequalities 

^ , I. ^ forallfc,/.e2;^ (C.26) 



(1 + |A; - /i|2)» " (1 + |yt - /i|2)r»-i 

71 



^ ^ ^ ^ max(l, -) for all rj > 0, z ^ . (C.27) 



^ V ' 

1 + 772; rj 
Step 1. For sdl k e Z'^, one has 

< AX^{k) ^ S%n{k) ; (C.28) 

here, as in (6.67), 

sx (k) ji+drxmr^' 

A being defined by (6.64). If A fulGUs (6.68), the above upper bound is such that 

SXn{k) = 0{jj^^) ^ for A; ^00. (C.29) 

The inequality AXn{k) > is obvious. To prove the rest we start from Eq. (C.26), 
implying 

AX^{k) ^ + j^^' ] Yl (l + |/,|2)n(l+|fc_/,|2)n-l 5 (^-30) 



setting 

we infer from (C.30) that 



li(k):^ inf \k-h\, (C.31) 

hez<i,\h\^A{k) 



1 + |A;|2 



where AS„ is defined following Eq. (C.12). To go on we claim that, for all k e Z'^, 

^^^)^{ A(A;)-I^l ifAS^I^I'. ^^-^^^ 

The above inequality is trivial if A(A;) < |A;|; if A{k) ^ |A;|, it follows noting that 

\h\ ^ A{k) implies \k - h\ ^ \h\ - \k\ ^ A{k) - \k\. 

Having proved (C.33), we insert it into (C.32); this gives the inequality 

AXn{k) ^ X{ky'-^ASn{Hk)) (C.34) 
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and substituting therein the upper bound (C.13) for AS„ we get the upper bound 
in (C.28). 

To go on, suppose A fulfills (6.68); then 
m > ^\k\ > \k\, X{k) < 1 ^ ;).|^|. ^ niax(l, ^— for \k\ ^ x (C.35) 
(the last inequality follows from (C.27), with z = |A;p and r] = {a — 1)^). So, 

^ 0{^) , X{k) = 0{l) forA;^oo; (C.36) 



A{k)-Vd \k\ 

inserting this in the definition of (5DC„, we infer Eq. (C.29). 
Step 2. With E„ as in (6.57), one has 

X'^{k) ^ E„ for A; ^ oo . (C.37) 

To prove this, we write 

X^k) = ""-kih) , (C.38) 

1 if ;2; e (0, +oo) , 
if ^ e (-00,0] . 



e{\k\/2-\ 


\h\)\k-h\ 


p(i + i^n— 1 




f|/iP)"(l + l 


\k-h\^Y 



Cnk{h):- (27r)^(l + |/iP)"(l ' I'-- ^l2^n ' ^(^) - 
For any fixed /i e Z*^, we have 



Cnfc(/^) ^ {27rY{l+\h\^r k^+oo; (C.39) 

this imphes 

V 7 T7T^ = , (C.40) 

if the limit /c — > oo can be exchanged with the sum over h. By the Lebesgue theorem 
of dominated convergence, the exchange is possible if Cnk{h) is bounded from above 
by a summable function of /i, uniformly in k] indeed this occurs, since 

^nk{h) ^ 1^ (TP])^ h,keZ<^. (C.41) 

Let us prove (C.41). The thesis is obvious if ^ 1^1/2, since in this case Cnd{k) = 0; 
hereafter we assume \h\ < \k\/2. First of all, we note that 
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c Uh) - \k-hni + \k\r-' 421 



2\n-l 



(27r)'^(l + + |A; - ' 

secondly, from \h\ < \k\/2 we infer \k — h\ ^ — \h\ ^ |A;|/2, whence 

1 / l + \k\^ \ 1 4^2^ 1 

CnfcW ^ ^27r)'^ Vl + \k\y4 ) (1 + \h\^r " (27r)'^ (1 + ^ ^ ^ 

(the last inequality follows from (C.27), with r] = 1/4 and z — 
Step 3. Suppose A fulfills (6.68); then, for A; — > oo, 

0{^) iin>d/2 + l, 



Kik) 



k\ 

og\k\. 



0{-^) iin = d/2 + l, ^0. (C.44) 

\k\'2n-d^ 



\k\^ 

OiTTT^) iid/2<n<d/2 + l 



To prove this, let us take any k E Z'^ such that \k\ ^ x- First of all, from (C.26) 
and A(A;) ^ P\k\ we infer 

DC"m<ii±^^ V I fC45) 

"^^^ {27rY , ^ (l+|/,|2)n(l + |^_/,|2)n-l ^'^'^^^ 

^ ^ ^eZd,|ik|/2^|/i|<A(ik) V ' I I ; V I I I ; 

(27r¥ ^ (l + |/i|2)"(l + |A;-/iP)"-i • 

For \h\ ^ |A;|/2 we have 1/(1 + \h\^) ^ 1/(1 + |fc|V4), whence 

^ " (27r)'^(l + |A;|V4)- ^ (i + |A; - /iRn-i ^^"^"^ 



l + |A;p 1 



(2.)^ VI + \k?N 1 + 1/^174 (1 + 1/^ - /^n^- 

4"-^ 1 1 
^ (2'nY 1 + |A;|2/4 ^ (1 + lA; - /iP)'^-! 

^ ^ ' ' ' heZd,|fc|/2^|fe|</3|fc| ^ ' ' ^ 

(the last passage uses again (C.27), with r\ = 1/4 and z = |A;p). Now, a change of 
variable h = q — k in the last sum gives 

^"^^^ ^ (27r)'^l + |A;|V4 ^ (1 + " ^^^'^^^ 

74 



To go on we note that, for all q e Z'^, 

\q-k\<(3\k\ \q\ ^ \q - k\ + \k\ < {(3 + l)\k\ . (C.48) 



Thus, 



§n-i{{P + mi) , 



1 + |A;|V4 

where S„_i is defined following Eq. (C.9). The A; — > oo behaviour of S„_i((/3+ 1)|A;|) 
is inferred from Eq. (C.ll); this function is multiplied by 4'*~^/(l + |A;p/4) = 

and the conclusion is Eq. (C.44). 
Step 4. Conclusion of the proof. For any cutoff A, the decomposition (C.21) and 
the bounds (C.28) give the inequalities (6.65) for /C„(A;), also implying the finiteness 
of this quantity for arbitrary k. 

Prom now on A is supposed to fulfill (6.68), and we consider the hmit A; — > oo. Then, 
the decomposition 0C„ = X'^ + and Steps 2, 3 give 



as claimed in (6.69); the other statement in (6.69), namely, 53C„(A';) = 0(l/|A';p"'~'^) 
0, is known from (C.29). The decomposition Kn = 3C„ + AIK„, and Eqs. (C.28) 
(6.69) also give 

lCn{k) ^ E„ . (C.50) 
Summing up, all statements in items (i)(ii) of this lemma are proved. □ 



D Appendix. Proof of Propositions 6.3 and 6.5. 

Proof of Proposition 6.3. We consider two vector fields v,w & H" with n > d/2, 

and proceed in two steps. 

Step 1. Proof of Eq. (6.71). We have 

v'' = ^h(^h , drw' = i irW^ee ; (D.l) 

this implies v»dw^ = v^drW^ — i Ylih te'Z'^^'^h^'r'^D^^h^^^ ™ vector form, 

v,dw = i Y {vh»t)w^{ehet) . (D.2) 
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On the other hand ChC^ — e/i+£/(27r)'^/^, so 

^•^^ = 72^Z1( (D.3) 



(27r)<^/2 

^ ^ kezd h,eezdh+e=k 



( [vh» {k - h)]wk-hj Cfc 



(27r)'^/2 



The term multiplying Sk in the last expression is the Fourier coefficient {v»dw)k; so, 
Eq. (6.71) is proved. 

Step 2. Proof that v.dw is in W'^ and fulGUs (6.73). For any k e Z*^, Eq. (6.71) 
implies 

^ E \vh\\k - h\\wk-h\ (D.4) 



Now, Holder's inequality | J2h (^hbh\'^ ^ I^/^P) (E/^ l^ftp) gives 

\{v.dw)k\'^ ^ CfcPfc , (D.5) 



|A;-/i| 

Cfe — 



2 



Pk := J] (1 + i/ir)i^/.r(i + - /^D-kfc-.r • 

hez<^ 

The last inequality imphes 

J](l+|^|2)"-i|(^.at/;)fe|2^ ^(l + |A;|2)"-ic,p, (D.6) 



On the other hand. 



(1 + \k\'Y-^Ck = /C„(A;) , sup (1 + \k\'Y-\n ^ K (D-7) 

with )Cn, Kn as in Eqs. (6.75) (6.74); the finiteness of ICn{k) for any k, and of its 
sup over k are known from Lemma 6.2. To go on, we note that 

E = ( E (1 + i^i'r i-'^i') ( E (1 + i^i'r i^'^i') - mi mi ■ m) 

feezrf hezd hezd- 
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Inserting Eqs. (D.7), (D.8) into (D.6), we see that Efcezd(l + \k\^T~^\{v.dw)k\'^ < 
+00, implying v»dw £ H""^. Furthermore, 

yielding Eq. (6.73). □ 

Proof of Proposition 6.5. This is a simple variant of the proof given for Propo- 
sition 6.3. One takes into account the following facts: if f,g E H",,, then their zero 
order Fourier coefficients are /o = 0, go = 0; furthermore, {f»dg)o = 0, since this 
function has mean zero by Lemma 6.4. 

□ 



E Appendix. Derivation of the NS equations (6.87). 

The NS equations in physical and adimensional units. 

In the space L of (oriented) lenghts we fix some positive lenght A, determining the 
size of the system in consideration. The "space domain" of the system is modelled as 
a torus L'^/ {27rXY; we write x for any point in this domain, and t for the "physical" 
time. The NS equations are 

p{v{t) + v{t), dv{t)) = -dp{t) + riAv{t) + k{x, t) (E.l) 

where v{t) : x v{x, t), p{t) : x h- >• p{x, t) are the velocity and the pressure fields, 
while k{t) : x i-^ k{x, t) is the density of external forces; in the above, ' indicates 
the partial derivative with respect to time t. The coefficients p > 0, rj > are the 
density and viscosity of the fluid, respectively. The velocity fleld v{t) is required to 
be divergence free, to fulfill the condition of incompressibility. 

One passes to the adimensional form introducing three functions i/(a;, t), 7r{x,t), 
K{x,t) (x G T"*, t e [0, T) C R) via the equations 

2 2 

v(x,t) — -^i/lxjt) , p(x,t) — -^^^n(x,t) , k(x,t) — -^^^K(x,t) (E.2) 
pA pA pA 

X rjt 

With these positions, Eq. (E.l) is equivalent to 

z>(t) + iy{t)»diy{t) = -dnit) + Az/(t) + K{t) . 
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Hereafter we formalise this adimensional version, specifying the necessary functional 
spaces. First of all, we suppose 

KeC^'\[0,+oo),ir-^) ; (E.3) 

secondly, we stipulate the following. 

E.l Definition. The incompressible NS Cauchy problem with initial datum vq G 

E[^+-^, in the pressure formulation, is the following. 

Find V e C([0, T), H^+i) n C^{[Q, T), M^-^), tt e C([0, T), H'') such that 

v{t) + v{t)* dv{t) = -dT:{t) + Ai/(i) + K{t) for t e [0, T) , i/(0) = vq (E.4) 
(for some T e (0, +00]). 

We note the following: the requirements u e C([0, T), M^+i) and tt e C([0, T), i/") 
are sufficient for the right hand side of the above differential equation to be in 

c([o,r),Hr^). 

The equivalence between the pressure formulation (E.4) and the Leray 
formulation (6.87) of the Cauchy problem. 

Let us introduce the function 

77 e C°'^([0, +00), H^') , t ^ r]{t) := £K{t) . (E.5) 

For convenience, we report here the Cauchy problem in the form (6.87): 

Find u G C([0,T),H^+i) nCi([0,T),H2-i), such that 

i>(t) = Au{t) - £(z/(t). dp{t)) + ri{t) for t e [0, T) , z/(0) = vq 
(for some T e (0, +00]). The above mentioned equivalence can be stated as follows. 

E.2 Proposition. A function v of domain [0, T) fulfills (6.87) if and only if there 
is a function tt with the same domain, such that (i^, tt) fulfills (E.4). 

Proof. Suppose a pair (i^, tt) fulfills (E.4), and apply the projector £ to both sides 
of the differential equation. We have £z/(t) = z^(t), and £ commutes with both the 
time derivative ' and the Laplacian A; finally, £,d7r{t) = 0. These facts yield Eq. 
(6.87). Conversely, suppose a function u fulfills (6.87) and define 

7 e C{[0, T),ir-^) , t ^ 7(i) := Au{t) - u{t). dv{t) + K{t) - v{t) ; (E.6) 

then from (6.87) one infers £7(t) = 0, i.e., 7(t) e W^-^ for all t. Let 

TT : t e [0, T) ^ 7r(i) := d'^-iit) (E.7) 

with as in (6.37); then tt G C([0,T)i/") (since maps continuously W^~^ into 
//"). One easily checks that (i/,7r) fulfills (E.4). □ 
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F Appendix. Proof of Eq. (6.100). 

Let us rephrase the definition (6.99) of ^ as 

m = m - m) , (f.i) 

having put 

C:[0,T)^UI, t^C{t) such thatCix, t) ■.= ri{x + h{t),t) . (F.2) 

Clearly, the thesis (6.100) follows if we prove that 

CeC''\[0,+oc),mi;') . (F.3) 

To this purpose, we consider the Fourier coefficients ^(t), Cfe(^) of ri{t), ({t) and 
note that (F.2) imphes 

(kit) = r]k{t) e^'=*'^(*) {k e t e [0, T)) . (F.4) 

Let t,t' e [0,T). We have 

Ck{t) - Ck{t') = ak{t, t') + h{t, t') (F.5) 

afc(t,0:=r7fcWe''*'^*')(e^'=-('^(*)-'^(*'»-l) , /3fc(t, e^'*'^*') (r7fc(t) - r/,(t')) ■ 
Therefore ||C(t) - C(i')lln-i = VEfcez^l + |fcP)"-'ICfc(t) - Cfc(t')P has the bound 

||C(i)-C(Olln-i^^(i, + ^(^,0 , (F-6) 



^(^'0 := /$](1 + .B{t,t') := /5](1 + |fc|2)-i|/3,(t,t')P . 

y feezd y feezd 

On the other hand, Eq. (F.5) and the elementary inequality |e*'^'^ — 1| ^ 1^1 bl (fo^ 
all y e R*^) give 

W{t,t')\ ^ \k\ |%(t)| \h{t) - h{t')\ , mt,t')\ = \r^,{t)-r^,{t')\ . (F.7) 

Inserting these bounds into the espressions (F.6) of yl(t, t'), B{t,t') (and using (1 + 

\k\^T-^\k\^ ^ (1 + lA;^)") we get 

\m - mwn-i ^ \m\\n m) - h{t')\ + im - v{t')\\n-i . (f.s) 

Now, let us consider any compact subset / of [0, +oo). Then, the assumptions (6.95) 
on 7] and the nature of h ensure the existence of constants Q, Mi, M2 such that 
\\v{t)\\n ^ Q, \h{t)-h{t')\ ^ Mi\t-t'\ and \\7j{t) -7]{t')\\n-i < M2\t-t'\ iort,t' G /. 
This imphes 

\\C{t)-C{t')\\n-i<{QMi + M2)\t-t'\ , (F.9) 

and (F.3) is proved. 
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G Appendix. Proof of Proposition 7.2, item (iv). 

Our aim is to derive Eq. (7.13); we write this as 

sup Af{t)^V2, (G.l) 
te[o,+oo) 



J\f : [0, +00) [0, +00) , t ^ 7V(t) := [ ds ii_{t - sje"" . 



(G.2) 



Eq. (G.l) will follow from Steps 1 and 2, giving separate estimates on M{t) for 
^ i ^ 1/4 and t > 1/4. 

Step 1. For ^ t ^ 1/4 one has N'it) < -\/2. In fact, with this range for t Eq. 
(7.10) for /i_ implies 

e2* /•* , e2* ds e^* ^ 1 r- 

Af{t) = -y= / ds , ^ -y= / . = -y= 2\/t ^ — < . 



Step 2. One has sup ^y^/^J\f{t) = In fact, for alH > 1/4, using again Eq. (7.10) 
for we get 

N{t) = ( / ds+ ds)ii_{t-s)e-' = V2(l-e^/^-*)+Ce-* = v^-(^e^/^-C)e-* , 
Jo it- 1/4 



-1/4 p3s' 

C := / ds' 



lo V2es' 

One has the estimate C ^ 0.6 < \/2e^^'^, implying M{t) < a/2. From the above 
expression for jV(i), we also get limt_^+oo ■A/'(^) = V^; these facts yield the thesis. □ 



H Appendix. The constants K2 and K4 in dimen- 
sion d = 3. 

The above constants are needed for the numerical examples in Section 10; the route 
to compute them is outlined in Lemmas 6.1, 6.2 and Proposition 6.3. 

Computing We can take for it any constant such that 

/sup ]C2{k) ^ K2 . (H.l) 

We have the bounds (6.65): X2{k) < IC2{k) ^ X2{k) + 6X2{k), with the explicit 
expressions (6.66) for X2{k) and (6.67) for SX2{k) (and = Zq therein). Both X2 
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and 5%2 are defined in terms of some cutoff function A2, that we choose in this way: 
A2(A;) := 24 if < 4, and A2(A;) := 6|A;| if |A;| ^ 4. 

The above setting can be employed to evaluate }C2{k) for some set of values of k; 
computations performed for all A;'s with \ki\ ^ 10 (i = 1, 2, 3) seem to indicate that 

sup }C2(k) = lim }C2(k) . (H.2) 
feezg 

On the other hand, according to (6.62), 

lim IC2{k) = E2 , (H.3) 

where E2 is the series (6.57) with n — 2 and Z'^ — Z^. To estimate this series, we 
use the bounds (6.58): §2(A) < E2 ^ S2(A) + 5§2(A) with S2(A) and S§2{^) as in 
Eqs. (6.59) (6.60); these depend on a cutoff A, to be chosen as large as possible to 
reach a good precision. Taking A = 250, we get 

0.03607 ^ E2 ^ 0.03934 , (H.4) 

which implies, taking square roots. 



0.1899 < /sup }C2{k) ^ 0.1984 . (H.5) 



Retaining only two meaningful digits, we take as a final upper bound for ^sup /C2 
the quantity 

K2 := 0.20 . (H.6) 

Computing K4. We can take for it any constant such that y^supj^ggd K,4{k) ^ K4. 

We use again the bounds (6.65) )C4{k) ^ X4{k) + SX4{k); the cutoff A4 defining X4 
and SX4 is chosen setting A^ik) := 10 if < 10/3 and A4(A;) := 3|A;| if \k\ ^ 10/3. 
Computing {X4 + 5X4) (A;) for \ki\ ^ 6 (i = 1,2,3) we obtain numerical evidence 
that supj^.g2g(3^4 + SX4){k) is attained at A; = (3, 0, 0). So, 

sup }C4(k) ^ (X4 + 5X4) (3, 0, 0) ^ 0.004383 . (H.7) 

We also have 

sup K:4{k) ^ 3^4(3, 0, 0) ^ 0.004382 , (H.8) 
and in conclusion, taking square roots. 



0.06619 ^ /sup IC4{k) ^ 0.06621 . (H.9) 
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Retaining only two meaningful digits, we take as a final upper bound for \/sup /C4 
the quantity 

K4 := 0.067 . (H.IO) 
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